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ABSTRACT 

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have 
determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, 
using the observed luminosity function and clustering of 40,696 0.2 < z < 1.0 red galaxies in Bootes. Half 
of ~ 10^^'^ h~^ Mq halos host a red central galaxy, and this fraction increases with increasing halo mass. 
We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, 
although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that 
the stellar mass contained within the red population has doubled since z — 1, with the stellar mass within 
red satellite galaxies tripling over this redshift range. In cluster mass halos (> 10^^ h~^ Mq) most of 
the stellar mass resides within satellite galaxies and the intra-clustcr light, with a minority of the stellar 
mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are 
proportional to halo mass to the power of ~ 0.35. We thus conclude that halo mergers do not always lead 
to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, 
the stellar masses of their central galaxies typically grow by only ~ 30%. 

Subject headings: galaxies: evolution - galaxies: luminosity function, mass function - galaxies: elliptical and 
lenticular, cD - (cosmology:) large-scale structure of universe 



1. INTRODUCTION 

In a universe where the bulk of the mass is in the 
form of coUisionless cold dark matter (CDM), galaxies 
reside within gravitationally bound halos of CDM par- 
ticles (White & Rees 1978). The motion of these parti- 
cles is effectively governed by gravity alone, so both the 
mass function and clustering of dark matter halos are 
predictable, albeit non-trivial, functions of redshift and 
cosmological parameters (e.g.. Press & Schechter 1974; 
Sheth & Tormen 1999; Jenkins et al. 2001; White 2002; 
Tinker et al. 2005). Gravitational collapse results in the 
growth of dark matter halos over cosmic time while dy- 
namical friction leads to the orbital decay and merging 
of substructure within halos. As a consequence, there is 
the expectation that the galaxies residing within these 
halos will also grow via mergers. 
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In a CDM cosmology, halos more massive than groups 
(10^'^ h~^ Mq) undergo rapid growth via merging be- 
tween z = 1 and the present day. If the stellar mass 
within every dark matter halo was always contained 
within a single central galaxy, we would expect the most 
massive galaxies to undergo rapid growth via mergers at 
z < 1 . Until relatively recently, the rapid growth of mas- 
sive galaxies at z < 1 was an almost universal prediction 
of CDM galaxy formation models (e.g.. White & Frenk 
1991; De Lucia et al. 2006). However, it is difficult to 
predict the rate at which substructure merges within ha- 
los (e.g., Taylor & Babul 2001; Benson et al. 2003; Taf- 
foni et al. 2003; Boylan-Kolchin, Ma, & Quataert 2008) 
and at least some simulations overestimate the rate of 
galaxy growth via mergers. 

As much of the stellar mass within groups and clus- 
ters resides within satellites (e.g., Sandage, Binggeli, & 
Tammann 1985; Lin & Mohr 2004) and the diffuse intra- 
cluster light (ICL; e.g., Kemp & Meaburn 1991; Gonzalez 
et al. 2000; Arnaboldi et al. 2002; Feldmeier et al. 2004; 
Zibetti et al. 2005), not all halo mergers funnel stellar 
mass directly into central galaxies. It is thus plausible 
that massive galaxies do not grow as rapidly as their host 
halos. The current generation of CDM galaxy formation 
models predict a range of assembly histories for massive 
galaxies, with rates of 2 < 1 stellar mass growth varying 
by up to a factor of 2 (e.g., Baugh et al. 2005; Bower 
et al. 2006; De Lucia et al. 2006; Naab et al. 2007). 
When plausible models produce widely varying predic- 
tions, there is a clear need for robust observations to test 



these predictions. 

There is compening observational evidence that galax- 
ies grow via mergers at z < 1. There are thousands 
of examples of merging galaxies in the literature, span- 
ning a broad range of galaxy types, mass and redshift 
(e.g., Arp & Madore 1987; Toomre & Toomre 1972; Con- 
selice et al. 2003; van Dokkum 2005; Rines, Finn, & 
Vikhlinin 2007; Mcintosh et al. 2007; Lotz et al. 2008). 
While there have been valiant attempts to measure the 
rate of galaxy growth using catalogs of merging galaxies 
or galaxy pairs, there are large uncertainties associated 
with the selection function of such objects and the time- 
scales for merging. At low redshift, measured rates of 
red galaxy stellar mass growth via mergers span from 
~ 2% per Gyr (e.g., Masjedi et al. 2006; Masjedi, Hogg, 
& Blanton 2007) to ~ 10% per Gyr (e.g., van Dokkum 
2005). Although the exact merger rate remains a matter 
of debate, clearly there is some assembly of red galaxies 
via merging at z < 1. 

The evolving space density of galaxies is a conceptu- 
ally simple method for measuring the growth of galaxy 
stellar masses. At the epoch where the most massive 
galaxies are assembled, their space density will increase 
with decreasing redshift. Although conceptually simple, 
in practice such measurements are difficult. Very lumi- 
nous galaxies are strongly clustered (e.g., Norberg et al. 
2002; Zehavi et al. 2005a), so uncertainties are domi- 
nated by the sample variance from large-scale structure 
rather than Poisson counting statistics. The most mas- 
sive galaxies typically possess strong 4000 A breaks and 
lack strong emission lines, which complicates both opti- 
cal photometry and spectroscopy of z > 0.8 galaxies. De- 
spite these difficulties, there is increasing evidence that 
the stellar masses of the most massive galaxies grow by 
50% or less at z < 1 (e.g., Bundy et al. 2006; Cimatti, 
Daddi, & Renzini 2006; Wake et al. 2006; Brown et al. 
2007; Scarlata et al. 2007). As a consequence, pure pas- 
sive evolution models without any mergers remain popu- 
lar (e.g., Cimatti et al. 2006; Wake et al. 2006; Scarlata 
et al. 2007), as they provide a better approximation of 
massive galaxy evolution than many simulations using 
CDM cosmologies. However, since there is compelling ev- 
idence for galaxy mergers, the rate of stellar mass growth 
must be non-zero. 

Measurements of the relationship between galaxy stel- 
lar mass and dark matter halo mass provide important 
clues as to how galaxies grow over cosmic time. For ex- 
ample, if galaxy stellar mass was directly proportional 
to host halo mass we would expect galaxies to grow as 
rapidly as their host dark matter halos during mergers. 
The relationship between galaxy stellar mass and halo 
mass has been explored, using various methods for esti- 
mating halo mass, including gravitational lensing (e.g., 
Hoekstra, Yee, & Gladders 2004; Mandelbaum et al. 
2006b), satellite galaxy velocities (e.g., Conroy et al. 
2007b), X-ray temperature (e.g., Lin & Mohr 2004), X- 
ray luminosity (e.g., Brough et al. 2008), group velocity 
dispersion (e.g., Brough et al. 2006), galaxy clustering 



(e.g., van den Bosch, Yang, & Mo 2003; Zheng, Coil, & 
Zehavi 2007) and matching the galaxy luminosity func- 
tion with the predicted halo mass function (e.g.. Vale 
& Ostriker 2004). For the most massive galaxies, these 
studies find that galaxy stellar mass scales as halo mass 
to the power of roughly a third. 

The halo occupation distribution (HOD; e.g.. Peacock 
& Smith 2000; Seljak 2000; Berhnd & Weinberg 2002; 
Cooray & Sheth 2002) describes the number and dis- 
tribution of galaxies within dark matter halos. As the 
clustering and space density of dark matter halos are pre- 
dictable functions of redshift, HOD models can be con- 
strained with the observed clustering and space density 
of galaxies. HOD models combine much of the relevant 
astrophysics (e.g., the halo mass function) with empiri- 
cal descriptions of how galaxies reside within dark matter 
halos. While such models cannot be considered complete 
descriptions of galaxy evolution, they do provide key in- 
sights into how galaxies grow over cosmic time. For ex- 
ample, if halo mergers efficiently funnel stellar mass into 
central galaxies, we would expect relatively little stellar 
mass to reside within satellite galaxies. 

In this paper, we determine the HOD using the mea- 
surements of the luminosity function and clustering of 
0.2 < z < 1.0 red galaxies selected from the Bootes 
field of the NOAO Deep Wide-Field Survey (NDWFS) 
and Spitzer IRAC Shallow Survey. This work builds 
upon our previous study of the HOD of the most massive 
(> 1.6L*) red galaxies in Bootes (White et al. 2007, see 
§5.3 for a summary). Red galaxies are an ideal popula- 
tion for measuring the HOD and testing the predictions 
of galaxy formation models. As red galaxies have low 
star-formation rates, the growth of individual galaxies at 
z < 1 should be dominated by merging (via dynamical 
friction) rather than star-formation, thus simplifying the 
comparison of CDM models and data. In CDM models 
the most massive halos undergo rapid growth via merg- 
ers at z < 1 , and if this results in rapid growth of central 
galaxies, then this should be observable within our sam- 
ple. The spectral energy distributions of red galaxies 
are dominated by old stellar populations, so the opti- 
cal luminosities of z < 1 red galaxies are tightly corre- 
lated with stellar mass. Although our halo masses rely 
upon CDM models of the space density and clustering 
of halos, we can compare our halo masses with those 
derived with other techniques (e.g.. X-ray temperature, 
weak lensing). In several galaxy formation models virial 
shock heating (e.g., Birnboim & Dekel 2003; Dekel & 
Birnboim 2006; Khochfar & Ostriker 2007) or feedback 
from Active Galactic Nuclei (AGNs; e.g.. Silk & Rees 
1998; Wyithe & Loeb 2003; Croton et al. 2006; Hopkins 
et al. 2006) truncates star-formation in halos above a 
critical mass, and this mass scale should be observable 
with a well constrained HOD. 

The structure of this paper is as follows. In 521 we de- 
scribe the NDWFS and IRAC Shallow Survey, the object 
catalogs and our red galaxy sample. We present measure- 
ments of the evolving luminosity function and luminosity 



density of red galaxies in fj3l In 21 we discuss the cor- 
relation function of red galaxies as a function of both 
luminosity and redshift, including comparisons with the 
literature. We discuss the HOD modeling of the clus- 
tering and space density of red galaxies in SJSJ including 
a discussion of the evolving relationship between galaxy 
luminosity and halo mass. In ^we introduce an illustra- 
tive analytic model of the red galaxy HOD, and discuss 
how stellar mass evolves within the central and satellite 
red galaxy populations. We summarize our key results 
in SjTl Appendix A provides an overview of the analytic 
method used to determine preliminary angular correla- 
tion function covariance matrices. Appendix B describes 
the mock galaxy catalogs, which we used to determine 
the uncertainties of the luminosity function and our final 
estimates of the correlation function covariance matrices. 
Throughout this paper we use Vega magnitudes and 
a flat cosmology with Um = 0.25, fib = 0.043, Hq = 
72 km s^^ Mpc^^, as = 0.8 and Us = 0.97. Our choice of 
cosmology is similar to the WMAP-3 cosmology (Spergel 
et al. 2007) and matches the 1 h~^ Gpc'^ dark matter 
simulation we used to generate mock galaxy catalogs. 
For the analytic calculation of the HOD (Q, we define 
halos as spherical objects where the mean density within 
the sphere is 200 times that of the background. For 
the mock catalogs, halos are identified using a Friends- 
of-Friends algorithm which is discussed in detail in Ap- 
pendix B. We denote base ten logarithms with log and 
natural logarithms with In. 

2. IMAGING, CATALOGS AND THE RED 
GALAXY SAMPLE 

We selected our red galaxy sample from imaging of 
the Bootes field by NDWFS and IRAC Shallow Sur- 
vey. Our sample is extremely similar to that of Brown 
et al. (2007), and we refer the reader to that paper for a 
more thorough description of the surveys, source detec- 
tion, photometry, photometric redshifts, rest-frame prop- 
erties, catalog completeness and sanity checks. The only 
significant difference between this work and Brown et al. 
(2007) is that we adopt a flat cosmology with il,„ = 0.25 
and Ho = 72km s^^ Mpc^^ rather than fl„i = 0.30 and 
Hq = 70km s^^ Mpc^^. This cosmology matches that of 
the 1 h~^ Gpc dark matter simulation we use to gener- 
ate mock galaxy catalogs (Appendix |B]) and is similar to 
the WMAP-3 cosmology (Spergel et al. 2007). 

2.1. IMAGING 

In this paper we utilize optical and infrared imaging 
of Bootes from the NDWFS (Jannuzi & Dey 1999) and 
Spitzer IRAC Shallow Survey (Eisenhardt et al. 2004). 
The NDWFS is an optical (BwRI) and near-infrared 
{K) imaging survey of two « 9.3 deg^ fields with tele- 
scopes of the National Optical Astronomy Observatory. 
We utilize the third NDWFS data releasq^ of optical 
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imaging with the MOSAIC-I camera on the Kitt Peak 
4-m telescope. To obtain accurate optical colors with 
fixed aperture photometry across the Bootes field, we 
have smoothed copies of the released images to a com- 
mon Point Spread Function (PSF) with a full width at 
half maximum of 1.35". The IRAC Shallow Survey is 
8.5 deg^ of contiguous imaging with the Infrared Array 
Camera (IRAC; Fazio et al. 2004) of Spitzer Space Tele- 
scope. The Bootes field was imaged at wavelengths of 
3.6, 4.5, 5.8, and 8.0 fim with a typical exposure time 
of 90 sec. We utilize the 3.6 and 4.5 /im imaging to re- 
move contaminants (e.g., stars) from our galaxy sample 
and (when combined with the NDWFS photometry) for 
photometric redshifts . 

2.2. SOURCE DETECTION AND PHOTOM- 
ETRY 

We detected sources using SExtractor 2.3.2 (Bertin & 
Arnouts 1996), run on the /-band images of the NDWFS 
third data release. For this paper, we only include ob- 
jects which are detected within nominal subfield bound- 
aries. The subfields have small overlaps and we remove 
duplicate object detections from these regions. To min- 
imize contamination of the catalog, we exclude regions 
surrounding very extended galaxies and saturated stars. 
Regions without good coverage from both the NDWFS 
and the IRAC Shallow Survey are also excluded. The 
final sample area is 6.96 deg^ over a 2.9° x 3.4° field-of- 
view. 

We measured aperture photometry for each object us- 
ing our own code. SExtractor segmentation maps were 
used to exclude flux associated with neighboring objects. 
We corrected the photometry for missing pixels (e.g., 
bad pixels) using the mean flux per pixel measured in 
a series of annuli surrounding each object. Uncertainties 
were determined by measuring photometry at ~ 100 po- 
sitions within 2' of the object position and finding the un- 
certainty which encompassed 68% of the measurements. 
The accuracy of the photometry was verified by adding 
artificial galaxies to our data, recovering them with SEx- 
tractor and measuring their photometry with our code 
(Brown et al. 2007). 

2.3. PHOTOMETRIC REDSHIFTS 

We determined redshifts for our galaxies using the em- 
pirical ANNz photometric redshift code (Firth, Lahav, & 
Somervillc 2003; CoUister & Lahav 2004). ANNz uses ar- 
tificial neural networks to determine the relationship be- 
tween measured galaxy properties and redshift. It does 
not use any prior assumptions about the shape of galaxy 
spectral energy distributions (SEDs), though it does as- 
sume the relationship between observed galaxy proper- 
ties and redshift is a relatively smooth function. 

The basis of our training set is galaxies with spectro- 
scopic redshifts in Bootes. The AGN and Galaxy Evo- 
lution Survey (AGES, C. S. Kochanek et al. in prepa- 
ration) has obtained spectroscopic redshifts of ~ 16000 
/ < 20 galaxies, while several hundred additional red- 



shifts are available from a variety of programs with 4 
and 8 m-class telescopes. At z > 0.6 and / > 20.5 there 
are relatively few galaxies with spectroscopic redshifts, 
so we added artificial galaxies to the training set by (at 
fixed redshift) interpolating in color and extrapolating in 
luminosity. We trained and measured photometric red- 
shifts using the 4" aperture photometry and the 2nd or- 
der moments of the /-band light distribution. Using 4319 
red galaxies with spectroscopic redshifts, we find the la 
uncertainties of our photometric redshifts are ~ 0.1 in 
redshift at / = 22, decreasing to ~ 0.03 at / = 19.5. 

2.4. REST-FRAME PROPERTIES 

To measure absolute magnitudes and rest-frame colors, 
we used maximum likelihood fits of Bruzual & Chariot 
(2003) SED models to the optical photometry. Through- 
out this paper we use solar metallicity models with a 
Salpeter (1955) initial mass function, a formation red- 
shift of z = 4, and a broad range of exponentially declin- 
ing star formation rates. As with all SED models and 
templates, these models do have errors but these have 
little impact upon our principal conclusions (Brown et al. 
2007, §6.5). 

The 4" aperture photometry captures 86% or less of 
the total flux. We corrected for the flux outside this 
aperture by assuming galaxies within our sample have 
de Vaucouleurs (1948) profiles truncated at 7 half-light 
radii. We adopt the same size-luminosity relation as 
Brown et al. (2007) and assume this relation undergoes 
luminosity evolution described by a Bruzual & Chariot 
(2003) T = 0.65 Gyr stellar synthesis model with a for- 
mation redshift of z — 4. While the size-luminosity re- 
lation of red galaxies has some scatter, this should have 
little impact upon our results as 4" aperture photometry 
captures most of the flux from z > 0.4 red galaxies (see 
Figure 5 of Brown et al. 2007). 

2.5. THE RED GALAXY SAMPLE 

The distribution of galaxy colors is bimodal (e.g., Hogg 
et al. 2004), and selection criteria for red galaxies typ- 
ically fall near the minimum between the red and blue 
galaxy populations (e.g., Madgwick et al. 2002; Bell et al. 
2004; Willmer et al. 2006). We apply a similar approach 
here, and use the following rest-frame color selection cri- 
terion: 

U -V>lAO- 0.25 - 0.08 x {My ~5logh + 20.0) 

-0.42 X (z - 0.05) + 0.07 x (z - 0.05)^ (1) 

Our criterion selects galaxies with rest-frame U—V colors 
within 0.25 magnitudes of the evolving color-magnitude 
relation of red galaxies. This criterion allows comparison 
with the recent literature and is similar to the criterion 
of Bell et al. (2004). We also apphed apparent color cuts 
to remove contaminants while not significantly reducing 
sample completeness, and these are discussed in detail in 
Brown et al. (2007). 

The completeness of the catalogs was verified by 
adding artificial galaxies with de Vaucouleurs (1948) 



profiles to copies of our data and recovering them with 
SExtractor. Our catalogs are more than 85% complete 
for / < 23.5 galaxies with half-light radii of 0.5" or less. 
We limit the absolute magnitude range in each of our 
redshift bins so we can determine accurate redshifts and 
have a sample completeness of 85% or more. Our final 
sample contains 40,696 red galaxies with photometric 
redshifts between z = 0.2 and z = 1.0. 

3. THE RED GALAXY LUMINOSITY FUNC- 
TION 

Our 0.2 < z < 1.0 red galaxy sample is very sim- 
ilar to that of Brown et al. (2007), except we adopt 
a slightly different flat cosmology with flm = 0.25 and 
Hq = 72 km s^^ Mpc^^. This increases the z = 1 lumi- 
nosity distance by 3% and the z < 1 comoving volume by 
10%. We select our sample using absolute magnitude cri- 
teria, so our sample size is 4% larger than that of Brown 
et al. (2007). While these changes are small, they are 
comparable to our random uncertainties, so we present 
revised measurements of the red galaxy luminosity func- 
tion and luminosity density here. 

Our methodology for measuring the luminosity func- 
tion is very similar to that of Brown et al. (2007). We 
measure the evolving luminosity function of red galaxies 
using both the binned 1/Knax estimator (Schmidt 1968) 
and maximum-likelihood fits of Schechter (1976) lumi- 
nosity functions. Our photometric redshift errors are 
relatively small and their impact upon the measured lu- 
minosity functions is estimated to be negligible (Brown 
et al. 2007). We split the sample into four photomet- 
ric redshift slices: 0.2 < z < 0.4, 0.4 < z < 0.6, 
0.6 < z < 0.8 and 0.8 < z < 1.0, which have comov- 
ing volumes of 8.1 x 10"^, 1.8 x lO^^, 2.9 x 10~3 and 
3.9 X 10~^ h~^ Gpc'^ respectively. The total comoving 
volume of 0.01 h~^Gpc^ is comparable to that of the 
z < 0.15 2dF Galaxy Redshift Survey (2dFGRS). 

Unlike Brown et al. (2007), we estimate the uncer- 
tainties for the luminosity function using multiple mock 
galaxy catalogs, whose construction is described in Ap- 
pendix [B1 We expect the uncertainties derived from 
mock catalogs to be more robust than those determined 
with subsamples of the data, as individual large-scale 
structures can span multiple subsamples. Fractional un- 
certainties derived from mocks are on the order of 10%, 
and can be 50% larger than those determined with sub- 
samples. 

To quantify the evolution of stellar mass within the red 
population, we measure the luminosity density (the lumi- 
nosity weighted integral of the luminosity function). To 
determine the i3-band luminosity density {jb) we sum 
over the galaxy catalog {J27=i ^0 ^^ ^^'^^ redshift bin, 
apply corrections for sample incompleteness, and divide 
by the volume of the relevant redshift slice. Galaxies 
fainter than our magnitude limits contribute up to '^ 15% 
of the luminosity density, and we include their contribu- 
tion using an analytic approximation of the HOD (^. 
Our measurements of the luminosity function and lumi- 



nosity density are summarized in Tables [T] and [2l 

The overall evolution of the luminosity function is very 
similar to that reported by Brown et al. (2007) and our 
principal conclusions remain unchanged. The luminosity 
density decreases by only 27 ± 20% between z = \ and 
= 0, while the i3-band luminosity of an aging stellar 
population fades by a factor of 3 over the same redshift 
range. We thus conclude that the stellar mass within 
the red galaxy population has increased by 118 ± 45% 
since z — \. Mergers of red galaxies (without accom- 
panying star-formation) do not increase the stellar mass 
contained within the red population. As there is little 
star-formation within the red population, the increase in 
stellar mass must result from stellar mass being trans- 
ferred from the blue galaxy population to the red galaxy 
population (e.g., Bell et al. 2004; Brown et al. 2007; 
Faber et al. 2007). 

The transfer of stellar mass from the blue galaxy pop- 
ulation into the red galaxy population will have little 
impact on the space density of the most luminous red 
galaxies unless significant merging takes place, as very 
massive blue galaxies are rare at z < 1 (e.g., Bell et al. 
2004). The evolving space density of luminous red galax- 
ies thus provides strong constraints on the rate of galaxy 
stellar mass growth via galaxy mergers. The bright end 
(Afs — Slogft, < —21) of the red galaxy luminosity func- 
tion fades by 0.9 i?-band magnitudes between z — 1 and 
z = 0. We would not expect to see any fading if mas- 
sive galaxies were being rapidly assembled via mergers 
at z < 1. However, the bright end of the luminosity 
function does not fade by the ~ 1.2 i?-band magnitudes 
predicted by passive stellar population synthesis models 
(e.g., Bruzual & Chariot 2003) and the luminosity evo- 
lution of the fundamental plane (e.g., van Dokkum & 
Stanford 2003; Treu et al. 2005). We thus conclude that 
there is ongoing assembly of massive galaxies at z < 1, 
albeit at a rate that only increases their stellar masses 
by ~ 30% between z — 1 and z = 0. 

4. MEASURING THE RED GALAXY COR- 
RELATION FUNCTION 

The spatial correlation function of galaxies is one of 
the principal constraints on the HOD. In particular, the 
clustering of galaxies on scales less than 1 h~^ Mpc is 
a strong function of the number of satellite galaxies re- 
siding within halos. We measured the clustering of red 
galaxies with the angular correlation function, and then 
used the Limber (1954) equation to compare models of 
the spatial clustering with our observations. Angular 
correlation functions were determined using subsamples 
of red galaxies selected as a function of photometric red- 
shift, luminosity and space density. A summary of the 
subsample properties, including galaxy counts and pho- 
tometric redshift ranges, is provided in Table [31 



4.1. THE ANGULAR CORRELATION FUNC- 
TION 

We measured the angular correlation function with the 
Landy & Szalay (1993) estimator: 



Coi9) = 



DD - 2DR + RR 
RR 



(2) 



where DD, DR, and RR are the number of galaxy- 
galaxy, galaxy-random and random-random pairs at an- 
gular separation 6 ± S9/2. The pair counts were deter- 
mined in bins spaced by 0.2 dex between 10" and 0.44°. 
For a field of finite size, estimators of the correlation 
function are subject to the integral constraint where 



i2)dnidn2 







(3) 



(Groth & Peebles 1977), where 612 is the angle separat- 
ing solid angle elements dQi and dfl2- The integral con- 
straint results in a systematic underestimate of the true 
clustering. If the number density fluctuations in the vol- 
ume are small, and the angular correlations are smaller 
than the variance within the volume, then to first or- 
der the correlation function is simply biased low by a 
constant equal to the fractional variance of the number 
counts. Often, to remove this bias, the term 
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(4) 



is added to Cj{9) where J7 is the survey area. The value 
of fi^tJn, where n is an estimate of the mean number of 
galaxies per unit area, is the contribution of clustering 
to the variance of the galaxy number counts (Groth & 
Peebles 1977; Efstathiou et al. 1991). 

Our final estimates of the integral constraint use the 
variance of the galaxy number counts measured from 
multiple mock galaxy catalogs (Appendix [B|. Before 
constructing the mock catalogs, preliminary estimates of 
the integral constraint were derived from power-law fits 
to the angular correlation function. The integral con- 
straint depends on the value of uj{9) on large angular 
scales. On such scales, the angular correlation function 
of mock catalog galaxies is constrained by the cluster- 
ing of dark matter halos while power-law extrapolations 
of the correlation function are strongly dependent on the 
measured value of the power-law index. Both mock cata- 
logs and power-law fits (where the power-law index of the 
spatial correlation function is ~ —2) provide comparable 
estimates of the integral constraint, with cjo — 0.02^(1'). 
The significance of the integral constraint reaches a max- 
imum of ~ 0.5cr for the largest angular scales of our high- 
est redshift bins. On the angular scales where most of our 
constraining power comes from, the integral constraint is 
negligible. 



4.2. THE COVARIANCE MATRIX 

To fit a functional form to the measured angular cor- 
relation function, we need a model covariance matrix to 



estimate uncertainties of our data points and the cor- 
relations between them. Commonly used methods for 
estimating the covariance matrix include mock catalogs, 
subsampling or resampling of the data, and analytic ap- 
proximations. 

For our preliminary measurements, we estimated co- 
variance matrices using the analytic approximations de- 
scribed in Appendix [Xj Analytic covariance matrices 
can be computed quickly and do not include random 
noise, but can have large systematic errors. The final 
results presented in this paper use covariance matrices 
determined with mock catalogs, whose construction we 
discuss in Appendix |B] The mock covariance matrices 
naturally include the effects of sample variance and shot- 
noise without making assumptions about the linearity of 
the underlying clustering or special relations of higher 
order functions. The mock catalogs are generated using 
a 1 h~^ Gpc^ dark matter simulation and an analytic 
approximation of how red galaxies reside within dark 
matter halos (^S]). The analytic approximation of the 
HOD was initially constrained with angular correlation 
functions determined with the analytic covariance matri- 
ces, and was then verified using mock catalog covariance 
matrices. 

4.3. THE GALAXY REDSHIFT DISTRIBU- 
TION 

The angular correlation function, uj{9), can be derived 
from the spatial correlation function ^(r, z), using the 
Limber (1954) equation: 



w{9)-- 




dN 



dz 



(5) 



where dN /dz is the spectroscopic redshift distribution, 
and r{6, z, z') is the comoving distance between two ob- 
jects at redshifts z and z' separated by angle on the 
sky. If the functional form of the spatial correlation 
function is known, the spatial correlation function can 
be determined with accurate measurements of the angu- 
lar clustering and a robust model of the galaxy redshift 
distribution. 

To model the redshift distributions of our subsamples, 
we use the Schechter luminosity functions listed in Ta- 
ble[2l Photometric redshift uncertainties result in objects 
being scattered into and out of our subsamples. To model 
this, we use the measured photometric redshift uncer- 
tainties from Table 4 of Brown et al. (2007), which are a 
function of both galaxy magnitude and redshift, and as- 
sume Gaussian random errors. We have not attempted 
to model gross photometric redshift errors as they are 
rare in our sample (see Figure 3 of Brown et al. 2007). 
If we did not include these uncertainties, we would un- 
derestimate the mean physical separation of galaxy pairs 
and underestimate £,(r) by approximately 20%. The pho- 
tometric redshift and model spectroscopic redshift distri- 
butions for four of our red galaxy samples are shown in 



Figure [T] As the spectroscopic redshift distribution of 
neighboring photometric redshift bins overlap, neighbor- 
ing bins are not completely independent of each other. 

Our model redshift distributions are a function of the 
measured photometric redshift uncertainties. As a con- 
sequence, if our photometric redshift uncertainties are in 
error our measurement of ^(r) will also be in error. To 
verify this is not the case, we checked that our estimate 
of ^ir) did not depend on the width of our redshift bins 
(Brown et al. 2003). When a redshift bin is very narrow, 
the width of the model redshift distribution is dominated 
by the photometric redshift uncertainties. Conversely, 
the model redshift distribution for a very wide redshift 
bin has relatively little dependence on the photometric 
redshift uncertainties. We found that our clustering es- 
timates did not vary when we modified the width of our 
redshift bins, consistent with there being no gross errors 
in our model redshift distributions. 

4.4. POWER-LAW MODELS OF THE COR- 
RELATION FUNCTION 

Power-laws provide a simple empirical approximation 
of the spatial correlation function, enabling quick com- 
parison with the literature. Also, unlike the HOD, power- 
law models of the correlation function do not depend on 
models of the space density and clustering of dark mat- 
ter halos. However, spatial correlation functions mea- 
sured with large redshift surveys do show significant de- 
partures from power-laws (e.g., Zehavi et al. 2004; Zheng 
et al. 2007). Also, power-laws are a purely empirical de- 
scription of the galaxy clustering and are not physically 
motivated. Power-law approximations of the correlation 
function are thus useful but should be treated with due 
caution. 

We approximate the spatial correlation function with 
the standard power-law parameterization: 



£,[r) ^ (r/roY 



(6) 



where r is the galaxy pair separation in comoving coor- 
dinates and ro is the spatial scale where ^(r) = 1. As 
a consequence, the angular correlation function is also 
parameterized by a power-law: 



w{9) =w(l') 



1-7 



(7) 



We use (j-'(l') rather than uj{l°) to parameterize the 
strength of the angular clustering, as the measured value 
aj(l°) strongly depends on the measured value 7. 

We plot examples of power-law fits to the angular cor- 
relation function in Figures [H and [31 along with HOD 
model fits (see fj5]) for comparison. While HOD models 
provide better fits to spatial correlation functions mea- 
sured with large redshift surveys (e.g., Zehavi et al. 2004; 
Zheng et al. 2007), we find comparable x^ values for 
HOD model and power-law fits to our data. This is not 
altogether surprising, as our sample is much smaller than 
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Fig. 1. — The photometric (left) and model (right) redshift distributions for the n brightest red galaxies selected down 
to a space density threshold of 10~^ h^ Mpc^"^. To model the galaxy redshift distribution, we used the luminosity 
functions from Tableland convolved them with the measured photometric redshift uncertainties. The mean redshift 
of the model galaxy redshift distribution is strongly constrained by the upper and lower limits of the photometric 
redshift bins and has a weak dependence on Schechter function parameters. If we did not account for photometric 
redshift errors when modeling the redshift distribution, we would overestimate the number of galaxy pairs with small 
physical separations and underestimate the spatial clustering of galaxies by ~ 20%. 
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Fig. 2. — The angular clustering of Mb — 51og/i < —19.5 red galaxies in Bootes as a function of redshift. Power-law 
fits to the data are shown with the grey line while the best-fit HOD models are shown with the black curves. The 
datapoints at large angular scales are highly correlated with each other, and contribute little to the power-law or HOD 
fits to the data. If we exclude the last three data-points when fitting HOD models to the data, the resulting HOD 
parameters change only marginally. 
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Fig. 3. — The angular clustering of red galaxies as a function of redshift, determined using samples of n brightest 
galaxies per unit volume down to a space density threshold of 10^"^ h^^ Mpc. Power-law fits to the data are shown 
with the grey line while the best-fit HOD models are shown with the black curves. 



comparable SDSS samples and (for a given sample) angu- 
lar correlation functions have lower signal-to-noise than 
spatial correlation functions. 

We determined the value of tq for each subsample using 
the measured values of w(l') and 7, models of the galaxy 
redshift distribution and the Limber (1954) equation. We 
provide a complete list of our power-law correlation func- 
tion parameters in Table [31 We best constrain the spatial 
correlation function on scales of ~ 1 h^^ Mpc, and on 
scales larger than ~ 5 ft.~^ Mpc our power-law models 
should be considered extrapolations. As we show in Fig- 
ure[31 the scale where power-law and HOD models of ^(r) 
equal unity often differ by ~ 1 /i~^ Mpc, so one should 
treat tq values from this work and the literature with 
caution. 

In Figure [5] we plot the spatial clustering oi z ^ 0.8 
red galaxies in Bootes (parameterized by tq) along with 
results from the recent literature at comparable red- 
shifts (Brown et al. 2003; Wilson 2003; Meneux et al. 
2006; Phleps et al. 2006; McCracken et al. 2007; Ross 
et al. 2007; Coil et al. 2008). Measurements from 
the CFHT UH8k survey (Wilson 2003) and the CFHT 
Legacy Survey (McCracken et al. 2007) are systemat- 
ically low, as their model redshift distributions do not 
account for photometric redshift errors. The clustering 
of red galaxies fainter than Mb — 5logh = —20 is roughly 
constant as a function of both luminosity and redshift, 
with rg ~ 5 /i~^Mpc. There is some evidence that the 
spatial clustering of luminous red galaxies is correlated 
with luminosity at high redshift, as one would expect if 
the most massive galaxies reside within the most mas- 
sive (and strongly clustered) dark matter halos. As halo 



masses are a product of HOD modeling, we will return 
to this point later in the paper. 

The clustering of red galaxies in Bootes, as a function 
of both luminosity and redshift, is plotted in Figure [6l 
For comparison we also plot low redshift measurements 
from the 2dFGRS (Norberg et al. 2002) and SDSS (Ze- 
havi et al. 2005b). We have used the model SO col- 
ors of Fukugita, Shimasaku, & Ichikawa (1995) to in- 
clude the 2dFGRS and SDSS measurements in Figure [SI 
Wc find that the clustering of red galaxies is clearly a 
strong function of luminosity for galaxies brighter than 
AIb — 5\ogh = — 20. This trend disappears or weakens for 
fainter galaxies. The faintest SDSS data-point suggests 
the correlation of tq with luminosity continues to fainter 
magnitudes, but the power-law approximation of this 
subsample has an unusually high value of 7 (2.46) which 
may explain the unusually low estimate of tq. The clus- 
tering of red galaxies as a function of luminosity clearly 
evolves with redshift. However, as we expect aging stellar 
populations to fade by 1.2 i3-band magnitudes between 
z = 1 and z = 0, much of this evolution is due to lumi- 
nosity evolution rather than the evolution of large-scale 
structure. 

5. HALO OCCUPATION DISTRIBUTION 
MODELING 

Halo occupation distribution (HOD; e.g.. Peacock & 
Smith 2000; Seljak 2000; Berlind & Weinberg 2002; 
Cooray & Sheth 2002) modeling is a powerful tool for 
understanding the clustering of galaxies and how galax- 
ies reside within dark matter halos. The principle of 
the HOD framework is to link galaxies to their host dark 
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Fig. 4. — Power-law and HOD models of the spatial correlation function, inferred from fits to the angular clustering of 
Mb — Blogh < —19.5 red galaxies at 0.4 < z < 0.6. The angular scale corresponding to a transverse comoving distance 
of r is shown at the top of the plot. On scales of '--^ 1 h~^Mpc, power-laws approximate HOD models of the correlation 
function. On larger scales, power-law models derived from our data should be considered extrapolations. Estimates 
of the spatial scale where £,{r) = 1 derived from HOD and power-law models can differ by 1 h~^ Mpc. 
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Fig. 5. — The spatial clustering of red galaxies at z ~ 0.8, including our results and those from the literature (Brown 
et al. 2003; Wilson 2003; Meneux et al. 2006; Phleps et al. 2006; McCracken et al. 2007; Ross et al. 2007; Coil 
et al. 2008). All the samples plotted here use red galaxies brighter than an absolute magnitude threshold, with 
the exception of NDWFS DRl and CFHTLS, which use absolute magnitude bins which are ±0.5 mag wide. Spatial 
clustering measurements from the CFHT UH8k survey and CFHT Legacy Survey are systematically low, as their 
model redshift distributions do not account for photometric redshift errors. The clustering of red galaxies fainter than 
Mb — 5logh = —20 is roughly constant with luminosity and redshift, with rg ~ 5 h~^Mpc. 
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Fig. 6. — The spatial clustering of red galaxies in Bootes as a function of absolute magnitude and redshift. For 
comparison, we also plot low redshift measurements from the 2dFGRS (Norberg et al. 2002) and SDSS (Zehavi 
et al. 2005b). All the samples plotted here use bins ±0.5 mag wide. The clustering of red galaxies brighter than 
Mb — 5logh = —20 is clearly a function of luminosity, and this trends weakens at fainter magnitudes. Compared to 
the other measurements shown here, the power-law fit to the faintest SDSS sample has an unusually high value of 7 
(2.46) and a correspondingly low estimate of tq. The clustering of red galaxies as a function of luminosity evolves, but 
this is largely due to the fading of aging stellar populations rather than evolution of large-scale structure. 



matter halos, whose formation and properties can be pre- 
dicted by both simulations and analytic methods. As the 
evolving space density and clustering of dark matter ha- 
los are predictable functions of redshift, one can deter- 
mine the galaxy HOD with the observed space density 
and clustering of galaxies. By measuring the HOD as a 
function of both galaxy luminosity and redshift, one ob- 
tains the evolving relationship between galaxy luminosity 
and host halo mass. One can thus determine how galaxy 
stellar masses are growing with respect to the masses of 
the dark matter halos in which the galaxies reside. 

5.1. HOD MODELING OF THE CORRELA- 
TION FUNCTION 

The conceptually simplest method of implementing 
HOD modeling is to populate halos identified in cosmo- 
logical A^-body simulations. We use this approach to 
generate mock galaxy catalogs, which we discuss in de- 
tail in Appendix [BJ The simulation provides halo masses 
and positions, while the HOD specifies the mean number 
and spatial distribution of galaxies within halos. 

The functional form of the HOD is typically motivated 
by the results of galaxy formation models, and it is phys- 
ically reasonable to split the HOD into central and satel- 
lite galaxy components (e.g., Kravtsovet al. 2004; Zheng 
et al. 2005). By definition central galaxies are found at 
the center of halos and there is only zero or one central 
galaxy per halo (although the mean number of centrals 



per halo can lie between zero and one). Multiple satel- 
lite galaxies can reside within a halo, and their spatial 
distribution is often assumed to follow the dark matter 
distribution. 

Once a set of HOD parameter values has been chosen, 
one can step through the halo catalog and populate each 
halo with galaxies. The HOD provides the probabilities 
that a halo will contain a central galaxy and a particular 
number of satellites. One can thus determine the number 
of galaxies in a given halo using these probabilities and 
a random number generator. The HOD also provides 
the probability of a satellite being at a particular radius 
from the halo center, so a random number generator can 
be used to assign satellite positions within halos. Once a 
catalog is populated with galaxies, one can mimic the ob- 
servations and statistics used to measure the properties 
of the galaxy population (in this case, the luminosity and 
angular correlation functions). One can repeat the pro- 
cess, exploring a range of HOD parameter values, until 
a best-fit model is found. 

While it is conceptually simple to populate a simulated 
halo catalog, it is computationally expensive. Analytic 
methods of HOD modeling have been developed which 
allow quick computation of the galaxy space density and 
clustering statistics (HOD; e.g., Peacock & Smith 2000; 
Seljak 2000; Berhnd & Weinberg 2002; Cooray & Sheth 
2002). These methods use descriptions of halo properties 
which have been calibrated against A^-body simulations. 
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We employ an analytic methodology which is extremely 
similar to that of White at al. (2007) and Zheng et al. 
(2007). We utilize analytic approximations for the halo 
mass function (Jenkins et al. 2001), the biased clustering 
of halos (Tinker et al. 2005), the profile of dark matter 
within halos (NFW; Navarro, Frenk, & White 1996), the 
concentration of halo profiles (Bullock et al. 2001) and 
the non-linear dark matter power-spectrum (Smith et al. 
2003). The galaxy space density is the integral of the 
halo mass function multiplied by the mean number of 
galaxies per halo. In the analytic calculation, the two- 
point correlation function is the sum of two terms, the 
one-halo term and the two-halo term. 

The one-halo term results from pairs of galaxies which 
reside within the same halo (intra- halo pairs). This term 
is dominant on small scales and, by definition, is sensitive 
to the fraction of galaxies which are satellites. Given an 
HOD model, the one-halo term is just the distribution of 
intra-halo galaxy pair separations as a function of halo 
mass convolved with the halo mass function. The calcu- 
lation of this term is usually decomposed into contribu- 
tions by central-satellite and satellite-satellite pairs. We 
assume satellite galaxies follow the dark matter distribu- 
tion within halos, which we model with NFW profiles. As 
a result, the spatial distribution of central-satellite pairs 
follows an NFW profile while that of satellite-satellite 
pairs follows an NFW profile convolved with itself. 

The two-halo term results from pairs of galaxies which 
reside within different halos, and this term dominates 
the correlation function on large scales. On the largest 
scales it is equal to the mass correlation function times 
the mean galaxy-weighted halo bias squared. On scales 
smaller than ^ 1 h~^ Mpc one must include prescriptions 
for halo exclusion (halos cannot reside within each other) 
and scale-dependent bias. In practice it is easier to calcu- 
late the power spectrum of the two-halo term rather than 
directly calculate the spatial correlation function, and we 
refer interested readers to Zheng (2004) and Tinker et al. 
(2005) for further details. 

5.2. HOD PARAMETERIZATION 

We use a five parameter model to describe the mean 
number of central and satellite galaxies (brighter than 
some luminosity) per halo. We plot an example of this 
model in Figure [7| and discuss the model in detail below. 
The functional form of this model is motivated by HODs 
observed in galaxy formation simulations (Berlind et al. 
2003; Kravtsov et al. 2004; Zheng et al. 2005). The 
mean number of central galaxies per halo is modeled with 



(iVee„(A/)) = I 



1 + erf 



logAf - logMn 



where erf is the error function 



erf (a;) = —= 



flogM 



dt 



(8) 



(9) 



At a halo mass of Mmin, 50% of halos host a central 
galaxy. If the relationship between galaxy luminosity and 



halo mass had no scatter, (iVccn) would be modeled by 
a step function. In reality this relation must have some 
scatter, resulting in a gradual transition from (Ncen) — 
to (Nccn) — 1, whose width we quantify with the param- 
eter CTlogA/- 

We approximate the mean number of satellite galaxies 
per halo with a power-law truncated at a threshold mass 
of Mq: 



(iV,at(Af)) = (Afccn(Af)) 



M -M, 



M' 



(10) 



The parameter M[ corresponds to the halo mass where 
{Nsiit(M)) ~ 1 when (as is the case here) M[ > Mq 
and M{ > Afmin- When a = 1 and Af > Mq, the mean 
number of satellites per halo is proportional to halo mass. 
The number of satellites in halos of a given mass is as- 
sumed to follow a Poisson distribution, which is consis- 
tent with theoretical predictions (Kravtsov et al. 2004; 
Zheng et al. 2005) and current observational constraints 
(Yang et al. 2005; Ho et al. 2007; Yang, Mo, & van den 
Bosch 2007). 

For each subsample of the red galaxy catalog we tested 
a range of plausible HOD models. For each HOD model 
we determined the galaxy space density and spatial corre- 
lation function using the analytic methodology described 
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Fig. 7. — The mean number of central and satellite galax- 
ies per halo, as defined by equations [H] and [TOl Plot- 
ted is the HOD of Mb — 51og/i < — 18 red galaxies at 
0.4 < z < 0.6. By definition halos host either zero or one 
central galaxy, with 50% of halos of mass Afmin hosting a 
central galaxy. The parameter txiogM quantifies the mass 
range where (A^ccn) transitions from ~ to ~ 1. Satel- 
lites can reside in halos more massive than Mq, and at 
mass of M[ the mean number of satellites per halo is 
~ 1. In very massive halos, the number of satellites is 
proportional to halo mass to the power of a. 
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above, and then determined the corresponding angular 
correlation function using the Limber (1954) equation. 
For each HOD model we estimated y^ values using the 
full covariance matrices and a space density prior includ- 
ing fractional uncertainties determined using mock cata- 
logs. 

To rapidly explore the plausible range of HOD param- 
eter space we applied the Markov Chain Monte Carlo 
method (MCMC; e.g., see Gilks, Richardson, & Spiegel- 
halter 1996). This method generates a list (or chain) 
of HOD parameters whose frequency in the chain traces 
the likelihood of that model fitting the data. It works 
by generating random HODs from a trial distribution 
and accepting or rejecting them based on the relative 
likelihood of the fit. We choose new models by perturb- 
ing the HOD parameters from the last accepted chain 
element by Gaussian offsets in the log of the relevant 
parameter. The step size and directions are determined 
from the covariance matrix of a previous run of the chain. 
We further restricted the HOD parameter space to those 
models with ciogM < 0.6. Each chain provides the HOD 
distribution and a set of models which provide good fits 
to the observations. 

We provide a summary of our best-fit HOD parame- 
ter values in Tables H] and [5l While Afmin and Af{ have 
small uncertainties, the other parameters are poorly con- 
strained. As the measured HOD parameter values are 
correlated with each other, Afmin and M{ as a function of 
luminosity would show less scatter if the other HOD pa- 
rameters had well determined values. For this reason, in 
Table [5] we provide fits of HOD models where only M^\^ 
and Af{ are free parameters while the other parameters 
are fixed at uiogM = 0.3, a = 1 and Mq — Mmin. The 
fixed parameters have values that are similar to those in 
Table H] and the predictions of simulations (e.g., Zheng 
et al. 2005; Seo, Eisenstein, & Zehavi 2007). 

5.3. THE RED GALAXY HOD 

Models of the correlation function derived from the 
HOD provide a robust estimate of the large-scale bias 
factor, hg (where ^g ~ ^^^dm), which we plot as a func- 
tion of absolute magnitude in Figure [51 HOD models 
better constrain the bias factor than power-law extrap- 
olations, which depend on the measured value of 7 and 
result in a varying bias factor on large scales. The bias 
factor increases with redshift, as the spatial clustering 
of red galaxies evolves slowly while the underlying dark 
matter distribution evolves rapidly. As noted previously, 
the bias (or clustering strength) is a strong function of 
luminosity for the brightest red galaxies while varying 
little for galaxies fainter than Mb — 51og/i = —20. 

In White et al. (2007), we discussed the HOD of space 
density selected samples of red galaxies, and we recap 
several key results of that study here. We selected the n 
most luminous red galaxies per unit volume down to a 
space density threshold of lO"'^" /I'^Mpc^^. If red galax- 
ies underwent pure luminosity evolution (PLE) without 
galaxy mergers, a space density selected sample would 



select the same (fading) galaxy population with redshift 
and the spatial clustering would increase with decreas- 
ing redshift due to gravitational collapse. As we show in 
Figure O we find little or no evolution of the spatial clus- 
tering of these galaxies with redshift, which is contrary to 
PLE. The evolving bias factor of these galaxies does not 
evolve in the same manner as analytic approximations 
(Fry 1996) or simulations (White et al. 2007) of PLE. 
The PLE simulation, normalized to the z = 0.9 obser- 
vations, overestimates the number of satellite galaxies at 
z = 0.5 by roughly 50%. A simple solution is to remove 
a third of the satellites, by having them undergo merging 
or disruption between z = 0.9 and z = 0.5 (White et al. 
2007). 

In Figure [10] we plot the values of Af,„in and M[ as a 
function of _B-band absolute magnitude. As one would 
expect, the most massive central and satellite galaxies 
can only reside in the most massive halos, so Mmin and 
M'l increase with luminosity. One can also see that the 
halo mass required to host a galaxy of fixed B-band lu- 
minosity decreases with increasing redshift. It is not ob- 
vious from Figure [TOj if this trend results from evolution 
of the relationship between galaxy stellar mass and halo 
mass, or the evolution of galaxy stellar populations. 
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Fig. 8. — The bias factor of red galaxies as a function 
of threshold absolute magnitude. Data points at a given 
redshift are correlated with each other as bright and faint 
galaxies can reside within the same large-scale structures 
and red galaxies selected by the brightest absolute mag- 
nitude thresholds are subsets of the red galaxies selected 
with fainter absolute magnitude thresholds. Red galaxies 
are highly biased tracers of the underlying dark matter 
distribution. As the spatial clustering of red galaxies 
does not evolve as rapidly as the underlying dark matter 
distribution, the bias increases with redshift. 
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Fig. 9. — The spatial clustering and bias factor of red galaxies as a function of redshift. The most luminous red galaxies 
were selected down to a space density threshold of 10"'^ h"^ Mpc""^. The value of tq was derived from power-law fits 
to the data while the bias factor was determined using HOD modeling of the clustering and space density of galaxies. 
If red galaxies underwent pure luminosity evolution without mergers, gravitational collapse would result in the spatial 
clustering increasing with decreasing redshift. A pure luminosity evolution model without mergers (Fry 1996), which 
is shown with the solid line, does not match our observations. 
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Fig. 10. — The relationship between red galaxy threshold absolute magnitude and host halo mass, parameterized with 
Mniin (left) and M[ (right). Afmin is the mass where 50% of halos host a central galaxy while M[ corresponds to the 
halo mass where the mean number of satellites is ~ 1. We determined Mmin and M{ using HOD modeling (with am, a 
and Mq fixed) of the observed clustering and space density of red galaxies in Bootes. Both Mmin and M[ increase with 
increasing luminosity, as one would expect if the most massive galaxies reside within the most massive dark matter 
halos. 
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Fig. 11. — HOD model parameters Afmin (left) and M[ (right) as a fmiction of threshold _B-band absolute magnitude 
plus 1.2z. By adding 1.2z to the absolute magnitudes, we compensate for the fading of red galaxy stellar populations, 
so the x-axes of our plots are proxies for stellar mass. Grey lines denote our analytic approximations of the HOD 
(given by Equations [12] and [T3|) . Both Afmin and M{ exhibit little or no evolution. We thus conclude that the 
relationship between red galaxy stellar mass and host halo mass undergoes little or no evolution at z < 1. As the 
masses of dark matter halos increase with decreasing redshift, we expect galaxy stellar masses to also increase with 
decreasing redshift. 



In Figured!] we compensate for the fading of red galaxy 
stellar populations by adding 1.2z to the i?-band abso- 
lute magnitudes. This correction is comparable to the 
observed evolution of the fundamental plane (e.g., van 
Dokkum & Stanford 2003; Treu et al. 2005) and is 
similar to the luminosity evolution of stellar populations 
models which reproduce the optical colors of red galax- 
ies (e.g., Beh et al. 2004; Brown et al. 2007). The value 
of Mb + 1.2z is effectively a proxy for stellar mass, al- 
though the exact conversion to stellar mass will depend 
on details of the stellar initial mass function and the star 
formation history. One can see in Figure [TT] that Afmin 
and M[ show little or no evolution as a function of red 
galaxy stellar mass. We thus conclude that the relation- 
ship between red galaxy stellar mass and host halo mass 
undergoes little or no evolution between z = 0.9 and 
z = 0.3. 

5.4. COMPARISON WITH HODs FROM THE 
LITERATURE 

An important consistency check is to compare our 
HODs with those from the literature. Unfortunately, 
such a comparison is not straightforward due to differ- 
ences in sample selection, adopted cosmology, and the 
assumed functional form of the HOD. To simplify the 
comparison, we plot HODs as a function of galaxy space 
density rather than attempting to model the equivalent 
absolute magnitude threshold in the f?-band. As space 
density is strongly correlated with luminosity, this pro- 
vides a simple and well defined way of comparing samples 
selected using a variety of techniques. 

HODs in the literature often use cosmologies or halo 



definitions that differ from the ones adopted here. To 
shift halo masses to our fim = 0.25 and as = 0.8 cosmol- 
ogy, we use the prescription of Zheng et al. (2002). For 
this comparison, we define halos to be spherical objects 
with a mean density 200 times that of the background, 
and we modify halo masses from the literature to be con- 
sistent with this definition. We compare characteristic 
quantities which are not sensitive to the details of the 
HOD parameterization. We use Afmin and Afi, where 
Afmin is the halo mass where (A^ccn) = 0.5 and Afi is the 
halo mass where (iVsat) = 1- 

In Figures [12] and [T^] we compare our HOD estimates 
of Afmin and Afi with those from the literature. The 
HODs from the literature are mostly derived from red 
galaxy samples, although we have included some HODs 
determined using samples including both red and blue 
galaxies (e.g., DEEP2, SDSS main). The various studies 
are in broad agreement, including measurements derived 
from the two-point correlation function (Zehavi et al. 
2005a; Phleps et al. 2006; van den Bosch et al. 2007; 
Zheng et al. 2007; Blake, CoUister, & Lahav 2008; Wake 
et al. 2008; Padmanabhan et al. 2008, Z. Zheng et al., in 
prep.), three-point correlation function (Kulkarni et al. 
2007), satellite galaxy velocities (Conroy et al. 2007b) 
and lensing (Mandelbaum et al. 2006a). The agreement 
between our measurements and those derived from lens- 
ing and satellite galaxy velocities is encouraging, as those 
measurements of halo masses do not depend on models 
of the halo mass function. 

The relationship between galaxy space density and 
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Fig. 12. — Mmin as a function of galaxy space density, where A/min is the halo mass where (-/Veen) — 0.5. Random 
uncertainties are on the order of 10%, with the exception of masses from satellite velocities, where uncertainties are 
~ 50%. Satellite halo masses have been revised downwards by 30% to compensate for spectroscopic incompleteness (as 
discussed by Conroy et al. 2007b). Our measurements are in broad agreement with the recent literature, which is to 
be expected as many of the Mmin measurements are constrained with models of the halo mass function. A power-law 
fit to the SDSS main sample (grey line) approximates Mmin as a function of both space density and redshift. This is to 
be expected, as the space density of ^ 10^^ h^^ Af© mass halos can be approximated by a slowly evolving power-law. 
Departures from this power-law are expected for high mass and high redshift halos. 
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Fig. 13. — M\ as a function of galaxy space density, where M\ is the halo mass where (Ngat) — 1- The grey line 
is a power-law fit to the HOD derived from the space density and 2-point correlation function of SDSS main sample 
galaxies. There is broad agreement between the various studies, although the scatter is larger than in Figure [T2l 
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A^min at z — can be approximated by a power-law 
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and this relationship appears to evolve slowly with red- 
shift. Again, this is to be expected, as the mass func- 
tion of dark matter halos less massive than 10^'^ h~^ Mq 
evolves slowly and can be approximated by nh{> M) oc 
M^^ . This relationship should not hold for the most 
massive galaxies at high redshift, since the top end of 
the halo mass function is not a power-law and rapidly 
evolves. 

While the results from the literature are in broad agree- 
ment, there is no consensus on how the relationship be- 
tween galaxy stellar mass and host halo mass evolves. 
Some studies find this relationship evolves (Yan, Madg- 
wick, & White 2003; Conroy et al. 2007b; Zheng et al. 
2007j3, while others find little or no evolution (Lin et al. 
2006; Phleps et al. 2006). Measuring HOD evolution 
is complicated by several factors, which may explain the 
varying conclusions in the literature. Redshift dependent 
selection effects can produce errors in the observed rate 
of HOD evolution. Studies of the evolving HOD may also 
be comparing (somewhat) different galaxy populations at 
different epochs. Finally, the relationship between blue 
galaxy stellar mass and host halo mass may evolve at 
z < \, while the same relation for > L* red galaxies 
undergoes little or no evolution. 

6. RED GALAXY EVOLUTION FROM HOD 
MODELING 

As we show in Figure [TTl there is little or no evolution 
of the relationship between red galaxy stellar mass and 
host halo mass at z < 1. Motivated by this, we intro- 
duce an analytic approximation of the HOD where the 
mean number of galaxies (of a given stellar mass) resid- 
ing within halos (of a given mass) does not evolve. The 
simplicity of this model allows one to easily discern im- 
portant relationships present in our data, including the 
correlation of central galaxy luminosity with host halo 
mass fi j6.2p and the contributions of central and satellite 
galaxies to the total stellar mass within a halo (i 36.3p . 
Our analytic approximation also predicts how the lumi- 
nosity function of central and satellite galaxies evolves 
over a broad redshift range (' ij6.4l and ij6.5p . 

Our approximation does not include all the relevant 
astrophysics (e.g., dynamical friction) and has empirical 
components (e.g., the functional form of (Nsat)), so it 
should not be confused with a complete model of galaxy 
evolution. We expect departures from this approxima- 
tion (particularly at high z) and these will be informa- 
tive. In this respect our approximation is somewhat anal- 



^Yan et al. (2003) find that the relationship between L* and host 
halo mass doesn't evolve, but L* does not correspond to a fixed 
stellar mass. Conroy et al. (2007b) find that the relationship 
between stellar and host halo mass does not evolve for all but the 
most massive galaxies. 



ogous to the illustrative and simple pure luminosity evo- 
lution models. That said, our model does include many 
of the key components required to describe the evolu- 
tion of z < 1 red galaxies, including the aging of stellar 
populations and the growth of dark matter halos. 

6.1. AN ANALYTIC APPROXIMATION OF 
THE HOD 

Our analytic approximation of the HOD is largely con- 
strained by the Mmin and M{ measurements presented 
in Figure 1111 These measurements span a limited lu- 
minosity range, so we use other measurements to con- 
strain the HOD of the most massive galaxies. The satel- 
lite fraction decreases with luminosity, with at least 80% 
of Mb — Blogh + 1.2z < —20 red galaxies being central 
galaxies. We thus constrain the HOD of the most lu- 
minous galaxies by assuming one very luminous galaxy 
per halo, assuming some scatter (uiogM = 0.3) in the re- 
lation between central galaxy luminosity and halo mass, 
and matching the cumulative luminosity function of these 
galaxies with the cumulative mass function of dark mat- 
ter halos. 

We approximate the HOD parameters My^in and M[ 
(for galaxies brighter than an absolute magnitude thresh- 
old Mb) with a series of power-laws, 

M„,i„(/i-iAf0) = lO"«5 

_|_2q11.95 ^ ^Q0A0x[-19~{MB-!ilogh+1.2z)] 
_|_2q13.70 ^ ^Ql.l5x[-21-{MB-!ilogh+1.2z)] 

(12) 



and 



Af{(/l-lM0) = 1012-70 X 100.11x[-17-(A/^-51og/.+1.2.)] 

i-[^q14.60 ^ -[^Q0.85x[-21-(A/i3-51og/i+1.2z)] 

(13) 

Each power-law component provides the relationship be- 
tween galaxy luminosity and halo mass for a particular 
luminosity range. For example, the last term of Equa- 
tion [12] shows halo mass is proportional to central galaxy 
luminosity to the power of 1.15 x 2.5 for the most lumi- 
nous galaxies. As we show in Figure llll this approx- 
imation reproduces the measured values of Mmin and 
A'I[ with an accuracy (RMS) of 8% and 7% respectively. 
As the other HOD parameters are poorly constrained 
by our observations, we fix their values to Mq = Mmin, 
ciogM ~ 0.3 and a = 1. 

We caution that Equations [12] and [13] are only ap- 
proximations of the true HOD. They also depend upon 
the assumed cosmology and the analytic approximations 
discussed earlier. It is also plausible that the HOD of red 
galaxies evolves very slowly at z < 1. That said, we do 
expect the broad brush-strokes of our model to be valid 
and it is these which we discuss below. 

Using the mean number of galaxies (brighter than a lu- 
minosity threshold) per halo and the halo mass function, 
one can determine the cumulative luminosity function 
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by estimating the space density of galaxies for a series 
of absolute magnitudes. Taking the derivative of this as 
a function of absolute magnitude, one obtains the con- 
ventional luminosity function. In Figure [Ml we compare 
the observed red galaxy luminosity function with our 
HOD model and Schechter function fits. The accuracy 
of the HOD model is comparable to that of Schechter 
functions which were individually fitted to the data in 
each redshift bin. The HOD model has a consistently 
higher space density of very luminous galaxies than the 
Schechter function fits, particularly at z < 0.6. While the 
faint-end of Schechter luminosity functions is a power-law 
with a constant index, the power-law index of the HOD 
model luminosity function varies at faint magnitudes. 

In Figure [15] we plot the red galaxy luminosity func- 
tion and the HOD model split by host halo mass. As 
one would expect, only the most massive halos can host 
the most massive galaxies. However, while some of the 
faintest red galaxies are central galaxies in lower mass 
halos (~ 10^^ h~^ Mq), a large fraction of low luminos- 
ity red galaxies reside (as satellites) in relatively massive 
halos. Galaxies in cluster mass halos (> 10^** h^^ Mq) 



represent a small fraction of the overall red galaxy pop- 
ulation. 

As we show in Figure [151 the shape of the red galaxy 
luminosity function varies with halo mass. In low mass 
halos, the shape of the galaxy luminosity function can be 
approximated by a Schechter function with a a > while 
in high mass halos a ~ — 1 . The value of a in high mass 
halos is very similar to the measured luminosity function 
of red cluster galaxies (e.g., De Propris et al. 2003). 

The luminosity function of group and cluster galaxies 
(> 10^'^ ft.^^ Mq mass halos) has a wiggle at ~ 2L*, cor- 
responding to the transition from a satellite dominated 
luminosity function to a central dominated luminosity 
function. Although predicted by theory (Zheng et al. 
2005) and observed in groups (Yang et al. 2007), the 
size and location of the wiggle depends on both uiogM 
and A/q, which are poorly constrained by our data. The 
luminosity function of red field galaxies is often approxi- 
mated by a Schechter function with a = —0.5, and in our 
HOD model this results from faint red galaxies residing 
in halos spanning a broad mass range. 
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Fig. 14. — Binned l/Vmax red galaxy luminosity func- 
tions and our HOD model of the luminosity function 
(solid lines). Random uncertainties for the 1/Vinax data 
points are on the order of 10%. A pure luminosity evolu- 
tion model, which fades by 1.2 i?-band magnitudes per 
unit redshift and is normalized to the z = 0.9 luminosity 
function, docs not match our observations. The accuracy 
of our HOD model is comparable to Schechter functions 
individually fitted to the data at each redshift. The con- 
tribution of central galaxies dominates the bright end of 
the luminosity function while the contribution of satel- 
lite galaxies increases with decreasing luminosity and de- 
creasing redshift. 
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Fig. 15. — Models of the luminosity function, split by 
halo mass, derived from our analytic approximation of 
the HOD. Clearly the most luminous red galaxies reside 
within the most massive halos. The least luminous red 
galaxies can be central galaxies in ~ 10"'^^ h~^ Mq halos 
or satellites in more massive halos. The luminosity func- 
tion of galaxies in high mass halos has a wiggle at ~ 2L* , 
which is predicted by theory (Zheng et al. 2005) and 
has been observed in galaxy groups (Yang et al. 2007). 
However, the amplitude and location of the wiggle are 
functions ciogM and Mq, which are poorly constrained 
by our observations. 
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6.2. THE GROWTH OF CENTRAL GALAX- 
IES 

In Figure [14] we plot the central and satellite compo- 
nents of the red galaxy luminosity function, determined 
using our analytic approximation of the HOD. Clearly 
the top end of the luminosity function is dominated by 
central galaxies. As we noted previously, halo mass scales 
as central galaxy luminosity to the power of 1.15 x 2.5, so 
the stellar masses of central galaxies scale as halo mass to 
the power of 0.35. We constrain the luminosity-mass re- 
lation for the brightest red galaxies using the observed lu- 
minosity function (^ and the predicted halo mass func- 
tion, so our luminosity-mass relation for these galaxies is 
similar to the z '--^ empirical model of Vale & Ostriker 
(2004, 2006) by construction. Our luminosity-mass rela- 
tion for the most massive red galaxies is broadly similar 
to others in the literature, where halo mass was deter- 
mined using X-ray temperature (Lin & Mohr 2004; Lin 
et al. 2006), X-ray luminosity (Brough et al. 2008), and 
weak lensing (Mandelbaum et al. 2006b). 

As the relationship between central galaxy stellar mass 
and halo mass is approximated by a non-evolving power- 
law with an index of 0.35, the most massive galaxies do 
not grow as rapidly as their host halos. While a merger 
of comparable mass halos presumably results in a merger 
of central galaxies, it is plausible that such mergers are 
not 100% efficient and dump stellar mass into the dif- 
fuse ICL (e.g., Gonzalez et al. 2000; Arnaboldi et al. 
2002; Feldmeier et al. 2004; Zibetti et al. 2005; Monaco 
et al. 2006; White et al. 2007; Conroy, Wechsler, & 
Kravtsov 2007c; Conroy, Ho, & White 2007a). Minor 
halo mergers may produce satellite galaxies instead of di- 
rectly funneling stellar mass into the central galaxy. As 
a result, while mergers may increase the mass of a halo 
by 100% at z < 1, the stellar mass of the central galaxy 
typically grows by only 30%. This is consistent with 
recent measurements of the luminosity function, where 
the most massive galaxies do not grow as rapidly as the 
most massive dark matter halos (e.g.. Brown et al. 2007). 
Semi-analytic models and simulations with rapid growth 
of massive galaxies may be overestimating the ability of 
halo mergers to funnel stellar mass into central galaxies 
(e.g., Taylor & Babul 2001; Benson et al. 2003; Taffoni 
et al. 2003; Boylan-Kolchin et al. 2008). 

Our approximation of the HOD and the observed lumi- 
nosity function of red galaxies (Figure I14p are inconsis- 
tent with pure passive evolution without mergers. How- 
ever, it remains plausible that the efficiency of central 
galaxy growth via halo mergers continues to decline with 
increasing halo mass (e.g., Cooray & Milosavljevic 2005). 
As the rate of major mergers with > 10^** h^^ Mq ha- 
los declines with decreasing redshift, the rate of massive 
galaxy growth will taper off. It is thus plausible that 
the space density of the most massive galaxies at low 
redshift can be approximated by pure passive evolution 
above some mass threshold (e.g., Bundy et al. 2006; 
Cimatti et al. 2006; Wake et al. 2006). 

Our analytic approximation of the HOD has a lower 



mass limit, so the cumulative luminosity function of cen- 
tral red galaxies only marginally increases with decreas- 
ing luminosity at L < 0.2L*. Halos far less massive than 
2Q11.85 ^-1 j^^^ presumably host blue central galaxies. 
The color bimodality of galaxies also shows evidence for 
a transition mass (e.g., Kauffmann et al. 2003), which 
is broadly consistent with our results. If star formation 
were truncated as function of galaxy stellar mass, one 
would not expect to observe red dwarf spheroidals. Thus, 
the bimodality of galaxy colors may be better described 
by a truncation of central galaxy star-formation in halos 
above a critical mass. Presumably there is a mechanism 
which prevents gas from cooling in high mass halos, thus 
suppressing central galaxy star-formation. The most 
plausible mechanisms in the current literature are virial 
shock heating (e.g., Birnboim & Dekel 2003; Khochfar 
& Ostriker 2007; Dekel & Birnboim 2006) and feedback 
from AGNs (e.g.. Silk & Rees 1998; Wyithe & Loeb 2003; 
Croton et al. 2006; Hopkins et al. 2006). 

6.3. STELLAR MASS IN CENTRAL AND 
SATELLITE GALAXIES 

In Figure [16] we plot the B-band luminosity con- 
tributed by red galaxies as a function of halo mass at 
z ^ 0.1, determined using our analytic approximation 
of the HOD. We also plot measurements of B-band lu- 
minosity versus mass, where the halo mass was deter- 
mined with weak lensing (Mandelbaum et al. 2006b), 
satellite galaxies (Conroy et al. 2007b) and X-ray tem- 
perature (Lin, Mohr, & Stanford 2004; Lin & Mohr 
2004). To include these measurements on our plot we 
assume Bvcga - r^B = 1-32 and i?vcga - -f'^Vega = 4. 
There is broad agreement between these measurements 
and our analytic approximation of the HOD, which we 
constrained with measurements of the space density and 
clustering of red galaxies. 

As we show in Figure [HI the fraction of S-band light 
contributed by satellite galaxies rapidly increases with 
halo mass. This is not unexpected, as the bulk of stellar 
mass within clusters does not reside within the brightest 
cluster galaxy (e.g., Sandage et al. 1985; Lin et al. 2004; 
Lin & Mohr 2004). Clearly the growth of these halos did 
not result in stellar mass being efficiently tunneled into 
the central galaxy. 

In high mass halos the total U-band luminosity is pro- 
portional to halo mass to the power of ~ 0.9 rather than 
1. This may indicate that mergers with very massive ha- 
los shift stellar mass out of the red galaxy population, 
perhaps into the ICL (e.g., Gonzalez et al. 2000; Arn- 
aboldi et al. 2002; Feldmeier et al. 2004; Zibetti et al. 
2005; Monaco et al. 2006; White et al. 2007; Conroy 
et al. 2007a, 2007c). 

6.4. THE EVOLVING LUMINOSITY DEN- 
SITY OF RED GALAXIES 

The luminosity weighted integral of the luminosity 
function, the luminosity density, is often used to esti- 
mate the evolution of the stellar mass within the red 
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galaxy population. In Figure[T7]we plot i3-band luminos- 
ity density measurements from Bootes and the literature 
(Madgwick et al. 2002; Bell et al. 2004; Blanton 2006; 
Willmer et al. 2006; Faber et al. 2007). Roughly 15% of 
the luminosity density is contributed by galaxies fainter 
than our magnitude limits, and we have applied correc- 
tions to the Bootes measurements using our HOD ap- 
proximation of the luminosity function. While there are 
some discrepancies between the various surveys, which 
we discuss below, there is a broad consensus that the B- 
band luminosity density of red galaxies evolves slowly at 
z < 1 (Bell et al. 2004; Brown et al. 2007; Faber et al. 
2007). 

In Figure [17] we over-plot the evolving luminosity den- 
sity derived from our analytic approximation of the 
HOD. As the stellar populations of red galaxies fade be- 
tween z = 1 and z = 0, a non-evolving luminosity density 
does not correspond to a fixed stellar mass density. To 
simplify the interpretation of the luminosity density, in 
the right-hand panel of Figure [IT] we divide the B-band 
luminosity density by io°'*^^'^^, so the y-axis is propor- 



I o 



o 

XI 

I 
m 



■ III 1 1 — 1 — 1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 1 1 ♦t: — 1 — X ly 1 1 1 

4*^ ■ 




/W- ..: 




i-j 




t/^T" 


' 




r ^ 




z~0.1 


/ 


HOD Total 


■/ 


HOD Centrals '■ 


■ ■/ 


L=M/(260h) 


/ 


¥ Sat. Velocities 


/ 


• Lansing 


A-i 


^ Cluster BCGs 


'1 1 1 1 1 . . 


♦ Clusters (Total) : 



10 



12 



10 



13 



10 



14 



10 



15 



Halo Mass (h"'' Mq) 



Fig. 16. — The B-band luminosity per halo contributed 
by all and central red galaxies at z — 0.1. Measure- 
ments of the relationship between galaxy luminosity and 
halo mass, derived from weak lensing (Mandelbaum et al. 
2006b), satellite galaxy velocities (Conroy et al. 2007b) 
and galaxy clusters (Lin et al. 2004; Lin & Mohr 2004), 
are in broad agreement with our analytic approximation 
of the HOD. (For clarity we do not show the uncertain- 
ties for individual clusters.) The B-band luminosity of 
all red galaxies is not directly proportional to halo mass 
(i.e., L — M/(260ft.)), which may indicate that stellar 
mass is lost from red galaxies during mergers with very 
massive halos (Monaco et al. 2006; Conroy et al. 2007c; 
White et al. 2007). 



tional to stellar mass. One can clearly see that the stellar 
mass contained within the red population doubles be- 
tween z = 1 and z = 0. As red galaxies have low rates of 
star-formation, presumably this mass is being transferred 
from the blue galaxies after a decline in star-formation 
(e.g., Beh et al. 2004; Brown et al. 2007; Faber et al. 
2007). 

In Figure [17] we show the contributions of central and 
satellite galaxies to the evolving luminosity density. If 
the stellar mass contained within the satellite and central 
populations remained fixed at their z = 1 values, the 
z = luminosity density of both populations would be 
~ 1.5 X 10^ h Lq Mpc~^ lower than what we derive from 
HOD modeling. The stellar mass contained within the 
red satellite population triples between z = 1 and z = 0. 
The bulk of z = red satellites were either central or 
blue galaxies within the past 7 Gyr. 

A naive interpretation of our results is that equal por- 
tions of new stellar mass are being transfered from the 
blue population to the red central and satellite popula- 
tions. Such an interpretation is risky, since satellites can 
merge with central galaxies, and central galaxies can be- 
come satellites in groups or clusters. One can imagine ex- 
treme scenarios where star-formation is only truncated in 
central or satellite galaxies. More sophisticated models, 
which track the evolution of individual halos and galax- 
ies over cosmic time will address these issues, and we will 
discuss such a model in an upcoming paper (M. White 
et al. in prep.). 



6.5. COMPARISON 
VEYS 



WITH OTHER SUR- 



In Figure [18] we plot measurements of the red galaxy 
luminosity function at low and high redshift (Madgwick 
et al. 2002; Blanton 2006; Bell et al. 2004; Willmer 
et al. 2006; Faber et al. 2007) along with the predic- 
tions of our HOD approximation and PLE. The PLE 
model, normalized to our z ~ 0.9 measurements, does 
not reproduce the red galaxy luminosity function. The 
PLE model underestimates the space density of ~ L* red 
galaxies and underestimates the luminosities of the most 
massive galaxies at low redshift. 

Our HOD model better approximates the galaxy lumi- 
nosity function than a PLE model. However, while our 
model reproduces the bright end of the luminosity func- 
tion, there are systematic offsets at low luminosities. As 
the most luminous red galaxies lie along the red sequence, 
their measured space density is insensitive to details of 
the galaxy selection criteria. In contrast, there is a sig- 
nificant population of < L* galaxies with colors falling 
between the red and blue galaxy populations. By shifting 
our rest-frame U — V selection criterion blueward by just 
0.1 magnitudes, the measured space density of L* red 
galaxies increases by 25% while the measured rg values 
decrease by 10%. Such a shift would remove the offset 
between our HOD approximation and the 2dFGRS and 
SDSS measurements. 

As the SDSS includes the Bootes field, we can directly 
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Fig. 17. — The i?-band luminosity density of red galaxies. The luminosity density values in the right panel have been 
divided by io0-*^i-2z^ go a fixed stellar mass corresponds to a horizontal line. Random uncertainties are on the order 
of 10% while systematic errors could be as high as 20% at z < 0.8. Pure luminosity evolution is ruled out by our 
observations and those from the literature. Our approximation of the HOD is shown with the black solid line, along 
with the contributions by central and satellite galaxies. The stellar mass contained within the satellite population 
rapidly increases with time. If the luminosity density of red galaxies rapidly declines at z > 0.8 (Bell et al. 2004; 
Willmer et al. 2006; Faber et al. 2007), the HOD must rapidly evolve. 
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Fig. 18. — Red galaxy luminosity functions for redshifts z ~ 0.1 and z ~ 1.1 (Madgwick et al. 2002; Bell et al. 2004; 
Blanton 2006; Willmer et al. 2006; Faber et al. 2007). For comparison, we also plot luminosity functions derived 
from our approximation of the HOD and a PLE model, which is normalized to our z = 0.9 data. For clarity we do 
not show the 2dFGRS and SDSS uncertainties, and we have increased the DEEP2 luminosities by 15% to correct for 
flux falling beyond the DEEP2 aperture (Brown et al. 2007). While our HOD model reproduces the bright end of the 
luminosity function, it fails for fainter galaxies where the measured space density is sensitive to details of the galaxy 
selection criteria. 
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compare our red galaxy sample with that of Blanton 
(2006). As we show in Figure [TH the Blanton (2006) 
red galaxy sample includes systematically bluer galaxies 
than our sample. This explains why Blanton (2006) mea- 
sures a higher space density of < L* red galaxies than 
Bootes. The rest- frame mq.i — .9o.i criterion of Blanton 
(2006) may be bluer than our U ~ V criterion, although 
other factors could be significant including the choice of 
galaxy SED templates and the method used to measure 
galaxy photometry. 

It is unclear if the HOD undergoes rapid evolution 
at z ~ 1. We find a significantly higher space den- 
sity of z ~ 0.9 red galaxies than either COMBO-17 or 
DEEP2. If the DEEP2 or COMBO-17 measurements of 
the z ~ 1 luminosity function are valid, the HOD must 
rapidly evolve at z ~ 1. As we noted previously, we do 
expect some evolution of the HOD as there is residual 
star-formation within the red galaxy population and the 
properties of dark matter halos of a given mass evolve. 
Why this evolution should rapidly accelerate at z ~ 1 is 
unclear. 

It is plausible that COMBO-17 and DEEP2 are un- 
derestimating the space density of z ~ 1 red galax- 




Redshift 

Fig. 19.— NDWFS Bw-R colors oi L > L* red galax- 
ies from the SDSS and this work. We use a 8" and 4" 
diameter apertures to measure the Byy — R colors of the 
SDSS and Bootes galaxies respectively, so the physical 
aperture size is comparable for both surveys. The SDSS 
0.05 < z < 0.15 red galaxy sample of Blanton (2006), 
selected using a rest-frame uo.i — go.i > 1.5 (AB) crite- 
rion, is systematically bluer than the evolving rest-frame 
U ~ V criterion used to select z > 0.2 red galaxies in 
Bootes. We thus expect the SDSS to measure a system- 
atically higher space density of ^ L* red galaxies than 
Bootes. 



ies. DEEP2 uses a non-evolving selection criterion for 
red galaxies, and this criterion must be redder than the 
evolving colors of an aging stellar population at high 
redshift. We thus expect DEEP2 to underestimate the 
space density of red galaxies beyond some (as yet un- 
determined) redshift. There is an ongoing collabora- 
tion between the Bootes and DEEP2 surveys to resolve 
this discrepancy, and we have obtained imaging of the 
DEEP2 Extended Groth Strip with the NDWFS filter- 
set to determine the impact of photometric errors on cur- 
rent luminosity function measurements. Unfortunately, 
COMBO-17 optical photometric redshifts have gross er- 
rors for galaxies beyond z = 1 (E. F. Bell, N. Taylor - 
private communication), resulting in a spuriously rapid 
decline in the space density of red galaxies. This sys- 
tematic error is being addressed with the addition of 
near-infrared photometry to the COMBO-17 photomet- 
ric redshifts. At present, it remains plausible that the 
relationship between galaxy stellar mass and host halo 
mass evolves slowly out to z ^ 1.5, and there is a clear 
need for robust space density and clustering measure- 
ments to test this conjecture. 

7. SUMMARY 

We have measured the past 7 Gyr of red galaxy growth 
within dark matter halos. To do this, we have measured 
the halo occupation distribution, which describes how 
galaxies reside within dark matter halos (as a function 
of halo mass). We have constrained HOD models using 
measurements of the evolving luminosity function and 
clustering of 0.2 < z < 1.0 red galaxies. Our sample of 
40,696 red galaxies, selected from 6.96 deg of imaging 
in Bootes, is far larger than comparable galaxy samples 
at these redshifts. 

We have measured the evolving luminosity function of 
red galaxies, and obtain results that are very similar to 
Brown et al. (2007). The bright end of the luminosity 
function fades with time, indicating that the most mas- 
sive galaxies are not rapidly growing via mergers. The 
bright end of the luminosity function does not fade as 
rapidly as an aging stellar population, and we conclude 
that the stellar masses of the most massive galaxies have 
grown by 30% over the past 7 Gyr. We find that the 
luminosity density of red galaxies evolves slowly with 
redshift, and the stellar mass contained within the red 
galaxy population has doubled between z = 1 and z = 0. 
As there is little star-formation within the red galaxy 
population, stellar mass must have been transfered from 
the blue galaxy population at z < 1. 

We have measured the angular correlation function of 
0.2 < z < 1.0 red galaxies, using samples selected as 
a function of photometric redshift, luminosity and space 
density. We evaluated power-law and HOD models of the 
spatial correlation function by using the Limber (1954) 
equation to determine corresponding angular correlation 
functions. To evaluate the quality of the models, we de- 
termined x^ values using the complete covariance ma- 
trices for the angular correlation function. Preliminary 
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covariance matrices were determined using analytic ap- 
proximations while our final results were determined us- 
ing covariance matrices derived from mock catalogs. Our 
mock catalogs reproduce the luminosity and correlation 
functions of red galaxies, and include large-scale struc- 
tures that are frequently observed in deep galaxy surveys. 

We find that the value of rg increases with luminos- 
ity for red galaxies brighter than Mb — 5\ogh — —20, 
while for fainter galaxies tq ~ 5 h^^ Mpc. While the 
spatial clustering of red galaxies as a function of lumi- 
nosity clearly evolves, this is largely due to the fading 
of aging stellar populations rather than the evolution of 
large-scale structures. 

We performed HOD modeling of the observed luminos- 
ity function and angular clustering of Bootes red galaxies. 
The large-scale bias factor increases with luminosity for 
red galaxies brighter than Mb — 51og/i = —20. We find 
that the bias of red galaxies increases with redshift, as the 
spatial clustering of galaxies does not evolve as rapidly as 
the underlying distribution of dark matter. The evolu- 
tion of the bias factor is inconsistent with pure luminosity 
evolution models without mergers (White et al. 2007), 
where the spatial clustering increases with time due to 
gravitational collapse. 

The relationship between red galaxy luminosity and 
host halo mass evolves, due to the fading of aging stellar 
populations. We do not observe any evolution of the 
relationship between red galaxy stellar mass and host 
halo mass with redshift. As halo masses grow at 2; < 1, 
we expect the stellar masses of central galaxies to also 
grow at z < 1. 

We find that the HOD of red central galaxies has a 
lower limit, with 50% of 10^^'^ h^^ Mq mass halos host- 
ing a red central galaxy. As one moves down further 
in halo mass, one presumably observes an increasing 
fraction of blue central galaxies. What truncates star- 
formation in central galaxies is unclear, but it appears 
to be strongly correlated with host halo mass. 

We can reproduce both the luminosity function and 
clustering of red galaxies in Bootes with an analytic ap- 
proximation of the HOD, which also simplifies the inter- 
pretation our of HOD measurements. In our model, the 
fraction of stellar mass within the satellite population in- 
creases with host halo mass. In cluster (> lO^** h'^ Mq) 
mass halos, more mass resides within satellites than cen- 
tral galaxies. We thus conclude that mergers with these 
halos do not always (efficiently) funnel stellar mass into 
central galaxies. The luminosities of central galaxies are 
proportional to halo mass to the power of ~ 0.35. As a 
result, if halo masses increase by 100% at 2; < 1, the stel- 
lar mass of the typical central galaxy grows by only 30% 
over the same redshift range. Massive galaxy assembly 
continues a,t z < 1, but the bulk of their growth took 
place at higher redshifts. 

The growth of stellar mass within the red population 
at z < 1 is reproduced by our analytic approximation of 
the HOD. The stellar mass contained within red satel- 
lite galaxies triples between z = I and z = 0, and they 



account for a third of the stellar mass within the red 
population by z ^ 0. Our analytic approximation of the 
HOD reproduces current measurements of the bright-end 
of the z > 1 red galaxy luminosity function, but there 
are offsets at lower luminosities. Either measurements 
of the 2 > 1 red galaxy luminosity function have large 
systematic errors, or the red galaxy HOD rapidly evolves 
at z > 1. 
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A. AN ANALYTIC APPROXIMATION FOR THE COVARIANCE MATRIX 

Our final estimates of the angular correlation function covariance matrices are derived from mock catalogs. However, 
one needs estimates of the HOD to generate such mocks, so preliminary fits of HOD models to the angular correlation 
function utilized analytic approximations of the covariance matrices. The basis of our analytic covariance matrices is 
the Gaussian approximation (e.g., Eisenstein & Zaldarriaga 2001): 

1 f°^ 
C^{6^,6,) = ^^ / KdK P^{K)Jo{Ke,)Jo{K0,) (Al) 

where Jq is a Bessel function and P2{K) is the (2D) angular power spectrum, 

/•oc 

P2{K)^2Tr 9de w{9)jQ{Ke) (A2) 

"'0 

This approximation is best suited to correlation functions where w(0) ^ 1 and underestimates the covariance of very 
strongly clustered objects. For a power-law power spectrum, uj{9) ex 9^"^ or P2{K) ex K°-~^ , the integral can be 
evaluated analytically in terms of = \j9i9j and r^ = 9il9j < 1: 

C.„,».*lSi:!l.^«-.I|±i) ,..,.- M-M,.') (A3, 

valid for 1 < a < 3/2. The power-series expansion of the (confluent) hypergeometric function, 2-Pl, converges rapidly. 
For cases where the power-law is (nearly) divergent we truncate the integral at an angular scale corresponding to a 
transverse comoving distance of 150/i~^Mpc. In practice, fits to our data are only marginally affected by the value 
of the angular correlation function on scales of more than 1°. If the angular correlation function bins have significant 
width. Equation lAll is modified to 

CU0^,O,) = I ,, \, 1 ( ,, \, 1 r%d0 r%'d9'CU0,0') (A4) 




(D. Eisenstein 2003, private communication) where 9i and 02 are the inner and outer radii of the bins. This can be 
rewritten as the single integral 

Cuj{9i,9j) 
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im2~&ii)^^' 

/ P^{K) [9^,2Jl{K9^,2) ' 9^.iJi{K9^.i)] [9j.2Ji{K9j.2) - 9j.iJi{K9,.i)] -—. (A5) 

Jo K 

The contribution of shot noise to the estimate of the covariance was included by adding the reciprocal of the sky 
surface density of galaxies (per steradian) to P2{K). However, the shot noise only dominates the covariance on scales 
of < 1' for the red galaxy sample. 

The analytic covariance matrix of Eisenstein & Zaldarriaga (2001) underestimates the uncertainties when lo{9) > 1. 
To mitigate this issue, we approximate the covariance matrix with, 

C'J9,,9,) = CUO^: 0j) + O.Olu;(0,)2-5 wherej > i (A6) 

where Cu, is determined using the method of Eisenstein & Zaldarriaga (2001) and uj{9j) is determined using the 
truncated power-law discussed earlier. The second term of this equation was determined by trial and error, and roughly 
approximates the near-diagonal elements of the Bootes covariance matrices derived from jack-knife subsamples and 
mock catalogs. While this term is a hack and should be treated with caution (and perhaps contempt), it does prevent 
highly correlated uj{9) measurements on small angular scales from skewing preliminary fits of power-laws and HODs 
to our data. 

In Figure[20]we plot the analytic and mock correlation matrices for one of our red galaxy samples. While the diagonal 
elements of the covariance matrices are comparable, the off-diagonal elements of the analytic covariance matrix are 
systematically less than those determined with mock galaxy catalogs. The uncertainties of power-law fits to the 
binned angular correlation function increase by up to 100% when we switch from analytic to mock covariance matrices. 
However, HOD parameters exhibit little change when we switch covariance matrices, as the HOD is constrained by 
the dark matter distribution, halo mass function and observed galaxy space density as well as the observed clustering 
of galaxies. 
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Angular Scale (Degrees) 

Fig. 20. — Analytic and mock catalog correlation matrices for the 0.6 < z < 0.8 Mb — 51og/i < —19.5 red galaxy 
sample. The correlation matrix is defined by C^^ {Oi ,0j)/ y^C^ {6i,9i) C^ {9j ,0j), where C^j is the covariance matrix of 
the angular correlation function. The near-diagonal elements of the two covariance matrices are in rough agreement, 
although the off diagonal elements show significant discrepancies. Despite this, our principal results do not change 
significantly if we switch from one set of covariance matrices to the other. 



B. MOCK CATALOGS 

We used mock galaxy catalogs to estimate the covariance matrix of the correlation function and the uncertainties 
for luminosity functions. We utilized a simulation which is very similar to that of White et al. (2007), but with a 
larger volume and lower mass resolution. The same cosmology is used as throughout this paper: Qm = 0.25 = 1 — VIa, 
r^b = 0.043, h = 0.72, Us = 0.97 and erg = 0.8. As in White et al. (2007), the hnear theory power spectrum for 
the initial conditions was computed by evolution of the coupled Einstein, fluid and Boltzmann equations using the 
code described in White & Scott (1996). Seljak et al. (2003) find that this code agrees well with CMBfast (Seljak & 
Zaldarriaga 1996). The simulation employed 1024^ particles of mass 10^^ h~^ Mq in a periodic cube of side 1 /i^^Gpc 
using a TreePM code (White 2002) with a Plummer-equivalent softening length of 35 h~^ kpc (comoving). A detailed 
comparison of this TreePM code with other codes can be found in Heitmann et al. (2007). 

We used a series of simulation outputs at 2 < 1 and for each output we generate a halo catalogs using the Friends- 
of-Friends (FoF) algorithm (Davis et al. 1985) with a linking length of 0.168 times the mean inter-particle spacing. 
This procedure partitions the particles into equivalence classes, by linking together all particle pairs separated by less 
than a distance b. The halos correspond roughly to particles with p > 3/(27r6^) ~ 100 times the background density. 
Our mass definition uses the sum of the particle masses in the halo, however to obtain better correspondence between 
our definition of halo mass and that implicitly defined by the mass functions of Sheth & Tormen (1999) and Jenkins 
et al. (2001) we rescaled the masses by M/Mfof = 1-1- 0.01 (InMfof — 23.5) where Mfof is the FoF mass in units of 
h~^MQ. With this redefinition the mass function in the simulation lies between those of Sheth & Tormen (1999) and 
Jenkins et al. (2001), differing from them by less than 10% in the mass range of interest. 

Mock catalogs were constructed by using an analytic approximation of the HOD (Sj6]) to assign central and satellite 
galaxies to halos. We determined the mean number of central red galaxies as a function of halo mass and used a 
random number generator to determine if a central galaxy was assigned to a particular halo. We determined the mean 
number of satellite galaxies as a function of halo mass, and used a random number generator and Poisson statistics to 
determine the number of satellites assigned to each halo. To determine the likelihood of a galaxy having a particular 
luminosity, we used conditional luminosity functions (number of galaxies per magnitude per halo) for central and 
satellite galaxies. These luminosity functions are a function of host halo mass and were derived from the cumulative 
luminosity functions provided by the analytic approximation of the HOD. Central galaxies were placed precisely at 
the center of dark matter halos while the distribution of satellite galaxies follows a spherical NFW profile. Random 
number generators were used to assign galaxy luminosities and satellite positions relative to halo centers. 

The 1 h~^ Gpc^ simulation was used to produce between 56 {z ^ 0.9) and 400(2; ~ 0.3) Bootes sized samples. The 
simulation cube was viewed parallel to one of the Cartesian axes, and each galaxy was assigned a spectroscopic redshift 
and corresponding apparent magnitude. Galaxies were then assigned photometric redshifts, assuming Gaussian random 
errors and the photometric redshift uncertainties reported by Brown et al. (2007), and the absolute magnitudes were 
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scattered accordingly. In Figure [211 we show the sky distribution of three 0.8 < z < 1.0 red galaxy mocks and the real 
0.8 < z < 1.0 red galaxy catalog. Individual large-scale structures at z ^ 0.9 are seen in both the Bootes field and our 
mock catalogs, and are relatively common in other deep surveys (e.g., Lubin at al. 2000; Nakata et al. 2005; Guzzo 
et al. 2007). While large structures of high redshift galaxies were once expected to be rare, this assumed rapid galaxy 
growth took place at z < 1 . 

Angular correlation functions were determined for each Bootes sized mock using the Landy & Szalay (1993) estima- 
tor, and these were used to derive covariance matrices. These covariance matrices include both sample variance and 
shot-noise without needing to assume the underlying fluctuations are Gaussian or higher-order moments follow some 
particular pattern. We also use the mocks to estimate the uncertainties for the galaxy luminosity functions. 

Are the mock catalogs consistent with the observed clustering of red galaxies in Bootes? To verify this we used the 
analytic approximation of the HOD to determine y^ values for the angular correlation functions measured from each 
mock. We then determined x^ in the same manner using the real Bootes field. In Figure [^ we plot the distribution 
of x^ values for the mocks and Bootes. With the possible exception of the faintest red galaxies in our highest redshift 
bin, the Bootes x^ values are comparable to those derived from the mocks. The analytic approximation of the HOD 
slightly underestimates the space density and overestimates the clustering of the faintest 0.8 < z < 1.0 red galaxies. 
This discrepancy may be caused by selection effects, as a small shift in the rest-frame selection criteria can alter the 
measured space density and clustering of < L* red galaxies by ~ 10%. As we discuss in ^QA[ such selection effects 
may contribute to the large scatter of z > 0.8 lunfinosity density measurements reported in the current literature. 
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Fig. 21. — The sky distribution of 0.8 < z < 1.0 red galaxies in Bootes (top-left) and three mock catalogs. Individual 
large-scale-structures are evident in both Bootes and the mocks, and are often found in deep galaxy surveys (e.g., 
Lubin et al. 2000; Nakata et al. 2005; Guzzo et al. 2007). While such high redshift structures were once thought to be 
rare, this assumed that massive galaxies and large-scale structure rapidly evolved at z < 1. Although one can imagine 
these structures grossly bias both luminosity function and clustering measurements, the number counts of the mock 
catalogs have a standard deviation of only 8%. 
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Fig. 22. — The distribution of x^ values for mock catalog correlation functions. The x^ values were determined using 
the angular correlation function measured from each mock and the analytic approximation of the HOD. Values of x^ 
for the Bootes field were determined in the same manner, and are shown with arrows. With the possible exception 
of the faintest red galaxies in our highest redshift bin, the Bootes x^ values are consistent with those derived from 
mocks. As we discuss in Appendix IbI the discrepancy between the 0.8 < z < 1.0 data and the mocks may be due to 
selection effects rather than an error in our approximation of the HOD. 
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Table 1 
Red Galaxy 1/Vmax Luminosity Function with mock catalog uncertainties. 

Absolute Luminosity Function {h Mpc^ mag~ ) 
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Table 2 
schechter function fits and the red galaxy luminosity density. 
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Table 3 
Red galaxy correlation function power-law fit parameters 
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5.6 ± 1.3 


0.37 


21.0 


< 


Ms 


- 51ogh 


< 


-20.0 





2 


< 


^ 


< 


0-4 


1071 


0.57 ± 0.14 


1.98 ± 0.19 


5.4 ± 1.5 


0.70 


19.0 


< 


Mp 


- 51ogh. 


< 


-18.0 





4 


< 


s 


< 


0.6 


3551 


0.46 ± 0.08 


1.93 ± 0.09 


6.5 ± 1.1 


1.22 


19.5 


< 


Mp 


- 51ogh. 


< 


-18.5 





4 


< 


z 


< 


0.6 


4270 


0.48 ± 0.07 


1.98 ± 0.08 


6.2 ± 0.9 


0.85 


20.0 


< 


Mp 


- 51ogh. 


< 


-19.0 





4 


< 


z 


< 


0.6 


4525 


0.47 ± 0.06 


2.00 ± 0.08 


5.9 ± 0.8 


0.66 


20.5 


< 


Mp 


- 51ogh. 


< 


-19.5 





4 


< 


z 


< 


0.6 


3875 


0.41 ± 0.06 


1.97 ± 0.10 


5.7 ± 1.0 


1.58 


21.0 


< 


Mp 


- 51ogh. 


< 


-20.0 





4 


< 


z 


< 


0.6 


2581 


0.44 ± 0.05 


1.95 ± 0.11 


6.0 ± 0.9 


0.67 


21.5 


< 


Mb 


- 51ogh. 


< 


-20.5 





4 


< 


^ 


< 


0.6 


1314 


0.64 ± 0.09 


2.02 ± 0.13 


6.7 ± 1.2 


0.96 


20.0 


< 


Mp 


- 51ogh 


< 


-19.0 





6 


< 


z 


< 


0.8 


5864 


0.24 ± 0.05 


1.91 ± 0.11 


5.2 ± 0.9 


0.24 


20.5 


< 


Mp 


- 51ogh 


< 


-19.5 





6 


< 


z 


< 


0.8 


5590 


0.25 ± 0.04 


2.08 ± 0.12 


4.6 ± 0.8 


1.04 


21.0 


< 


Mp 


- 51ogh 


< 


-20.0 





6 


< 


z 


< 


0.8 


4158 


0.28 ± 0.04 


2.11 ± 0.10 


4.7 ± 0.6 


0.55 


21.5 


< 


M^ 


- 51ogh 


< 


-20.5 





6 


< 


2 


< 


0-8 


2421 


0.42 ± 0.05 


2.17 ± 0.12 


5.3 ± 0.8 


1.09 


22.0 


< 


Mp 


- 51ogh. 


< 


-21.0 





6 


< 


z 


< 


0.8 


1090 


0.54 ± 0.11 


2.06 ± 0.20 


6.4 ± 1.5 


0.40 


20.5 


< 


Mp 


- 51ogh. 


< 


-19.5 





S 


< 


z 


< 


1.0 


7654 


0.19 ± 0.03 


1.93 ± 0.07 


4.9 ± 0.5 


0.36 


21.0 


< 


Mp 


- 51ogh. 


< 


-20.0 





S 


< 


z 


< 


1.0 


6534 


0.21 ± 0.03 


2.00 ± 0.10 


4.9 ± 0.6 


0.84 


21.5 


< 


Mp 


- 51og;i 


< 


-20.5 





8 


< 


z 


< 


1.0 


4147 


0.23 ± 0.04 


1.94 ± O.IS 


5.2 ± 1.0 


1.14 


22.0 


< 


Mp 


- 51ogh 


< 


-21.0 





8 


< 


z 


< 


1-0 


1992 


0.43 ± 0.06 


1.84 ± 0.14 


7.8 ± 1.1 


1.40 



Wide lun 


linosity 


ran 


go 


Mb 


— 51og/i 


< 


-17.5 


Ms 


- 51ogh 


< 


-18.0 


Mb 


- Slogh 


< 


-18.5 


Mb 


- 51ogii 


< 


-19.0 


Mb 


- Slog), 


< 


-19.5 


Mb 


- Slogh 


< 


-20.0 


Mb 


- Slogh 


< 


-18.0 


Mb 


- 51ogii 


< 


-18.5 


Mb 


- 51ogh 


< 


-19.0 


Mb 


- 51ogii 


< 


-19.5 


Mb 


- Slogh 


< 


-20.0 


Mb 


- Slogh 


< 


-20.5 


Mb 


— 51og?i 


< 


-19.0 


Mb 


- Slogh 


< 


-19.5 


Mb 


- Slogh 


< 


-20.0 


Mb 


- Slogh 


< 


-20.5 


Mb 


- Slogh 


< 


-21.0 


Mb 


— Slogh 


< 


-19.5 


Mb 


- Slogh 


< 


-20.0 


Mb 


- Slogh 


< 


-20.5 


Mb 


- Slogh 


< 


-21.0 



0.2 < 


z 


< 


0-4 


6027 


0.66 ± 0.11 


1.97 ± 0.10 


6.2 ± 1.2 


0.2 < 


z 


< 


0.4 


5275 


0.67 ± 0.10 


1.96 ± 0.09 


6.3 ± 1.1 


0.2 < 


z 


< 


0.4 


4419 


0.63 ± 0.11 


1.98 ± 0.11 


5.9 ± 1.2 


0.2 < 


z 


< 


0.4 


3385 


0.59 ± 0.11 


1.96 ± 0.12 


5.8 ± 1.2 


0.2 < 


z 


< 


0.4 


2189 


0.65 ± 0.13 


2.01 ± 0.13 


5.7 ± 1.3 


0.2 < 


^ 


< 


0.4 


1157 


0.68 ± 0.17 


2.03 ± 0.18 


5.6 ± 1.5 


0.4 < 


3. 


< 


0.6 


11117 


0.50 ± 0.05 


1.99 ± 0.06 


6.2 ± 0.7 


0.4 < 


z 


< 


0.6 


9541 


0.51 ± 0.06 


2.01 ± 0.06 


6.1 ± 0.7 


0.4 < 


z 


< 


0.6 


7566 


0.51 ± 0.05 


2.04 ± 0.06 


5.9 ± 0.7 


0.4 < 


z 


< 


0.6 


5271 


0.51 ± 0.06 


2.04 ± 0.08 


5.9 ± 0.8 


0.4 < 


z 


< 


0.6 


3041 


0.55 ± 0.07 


2.04 ± 0.10 


6.1 ± 0.9 


0.4 < 


^ 


< 


0.6 


1396 


0.68 ± 0.10 


2.16 ± 0.14 


5.9 ± 1.1 


0.6 < 


z 


< 


0.8 


11154 


0.28 ± 0.04 


2.07 ± 0.09 


4.9 ± 0.7 


0.6 < 


z 


< 


0.8 


8321 


0.31 ± 0.04 


2.09 ± 0.08 


5.0 ± 0.6 


0.6 < 


z 


< 


0.8 


5290 


0.33 ± 0.05 


2.17 ± 0.10 


4.8 ± 0.6 


0.6 < 


z 


< 


0.8 


2731 


0.41 ± 0.05 


2.21 ± 0.13 


5.1 ± 0.8 


0.6 < 


^ 


< 


0.8 


1132 


0.57 ± 0.12 


2.09 ± 0.20 


6.4 ± 1.5 


0.8 < 


g 


< 


1.0 


12398 


0.24 ± 0.03 


2.03 ± 0.07 


5.1 ± 0.5 


0.8 < 


z 


< 


1.0 


8619 


0.24 ± 0.03 


2.08 ± 0.08 


4.9 ± 0.5 


0.8 < 




< 


1.0 


4744 


0.27 ± 0.04 


2.08 ± 0.13 


5.1 ± 0.8 


0.8 < 


z 


< 


1.0 


2085 


0.44 ± 0.06 


1.87 ± 0.14 


7.7 ± 1.1 



0.86 
0.87 
0.84 
0.82 
0.68 
0.87 

1.29 
1.11 
1.11 
1.02 
0.81 
1.10 

0.73 
0.80 
0.45 
0.74 
0.57 

0.81 
0.79 
1.26 
1.27 



Space density selected sample 



-3.0^3 
-3.0 
-3.0 
-3.0 



Mpc~ 

Mpc~ 

Mpc~ 

h-^Mpc" 



-3.5,3 



^h^Mpc" 



0.2 < s < 0-4 


770 


0.78 ± 0.20 


2.05 ± 0.20 


5.8 ± 1.6 


0.99 


0.4 < s < 0.6 


1836 


0.57 ± 0.09 


2.08 ± 0.13 


5.9 ± 1.0 


1.44 


0.6 < s < 0.8 


2813 


0.42 ± 0.05 


2.14 ± 0.13 


5.4 ± 0.8 


0.88 


0.8 < s < 1.0 


3802 


0.34 ± 0.04 


2.03 ± 0.14 


6.0 ± 0.9 


1.56 


0.6 < s < 0.8 


955 


0.56 ± 0.14 


2.12 ± 0.27 


6.2 ± 1.9 


0.75 


0.8 < z < 1.0 


1233 


0.48 ± 0.10 


1.98 ± 0.18 


7.3 ± 1.4 


0.71 
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Table 4 
hod models of the luminosity and correlation functions of red galaxies. 



CO 
CO 



Photomctri. 
Rcdshift 



Spa 



X-'/d.o.f- 



'b 


— Slogh 


< - 


-17, 


5 





2 


< 


z 


< 


0.4 


6.98 ± 0.92 


X 


10' 


-a 


11.96 ± 0.08 


0.27 ± 0. 


19 


8.57 ± 2.76 


12.78 ± 0.13 


1.01 ± 0.13 


1.41 ± 0. 


.06 


0.77 


' B 


- Slogh 


< - 


-18. 








,2 


< 


z 


< 


0.4 


6.11 ± 0.80 


X 


10" 


-3 


12.00 ± 0.08 


0.30 ± 0. 


18 


10.60 ± 2.89 


12.84 ± 0.16 


0.99 ± 0.15 


1.40 ± 0. 


,06 


0.85 


'b 


- 51ogh 


< - 


-18. 


S 





.2 


< 


z 


< 


0.4 


5.14 ± 0.69 


X 


10" 


-3 


12.05 ± 0.08 


0.33 ± 0. 


18 


6.91 ± 4.39 


12.98 ± 0.14 


1.00 ± 0.14 


1.38 ± 0. 


,06 


0.93 


' B 


- Slogh 


< - 


-19. 








2 


< 


z 


< 


0.4 


3.96 ± 0.53 


X 


10" 


-3 


12.11 ± 0.08 


0.31 ± 0. 


19 


9.38 ± 3.62 


13.21 ± 0.13 


1.06 ± 0.19 


1.38 ± 0. 


.07 


1.02 


' B 


- Slogh 


< - 


-19. 


s 





2 


< 


z 


< 


0.4 


2.59 ± 0.36 


X 


10" 


-3 


12.27 ± 0.08 


0.28 ± 0. 


19 


10.73 ± 2.43 


13.46 ± 0.18 


1.03 ± 0.27 


1.40 ± 0. 


.08 


0.96 


'b 


— Slogh 


< - 


-20 








2 


< 


' 


< 


0.4 


1.42 ± 0.20 


X 


10" 


-3 


12.50 ± 0.09 


0.30 ± 0. 


19 


9.36 ± 2.64 


13.82 ± 0.18 


1.15 ± 0.33 


1.47 ± 0. 


.07 


1.17 


' B 


- Slogh 


< - 


-IS. 








4 


< 


z 


< 


0.6 


5.75 ± 0.53 


X 


10" 


-3 


12.02 ± 0.07 


0.35 ± 0. 


17 


11.58 ± 0.94 


12.81 ± 0.12 


1.07 ± 0.10 


1.57 ± 0. 


.03 


0.56 


' B 


- Slogh 


< - 


-18 


s 





4 


< 


z 


< 


0.6 


4.92 ± 0.46 


X 


10" 


-3 


12.06 ± 0.07 


0.32 ± 0. 


19 


11.87 ± 1.26 


12.86 ± 0.25 


1.06 ± 0.19 


1.57 ± 0. 


.03 


0.58 


'b 


- Slogh 


< - 


-19. 








4 


< 


z 


< 


0.6 


3.89 ± 0.37 


X 


10" 


-3 


12.14 ± 0.07 


0.32 ± 0. 


18 


12.38 ± 0.63 


12.86 ± 0.28 


0.94 ± 0.21 


1.59 ± 0. 


.03 


0.36 


'b 


- Slogh 


< - 


-19. 


s 





4 


< 


z 


< 


0.6 


2.71 ± 0.27 


X 


10" 


-3 


12.29 ± 0.07 


0.37 ± 0. 


18 


10.59 ± 3.06 


13.23 ± 0.15 


1.10 ± 0.16 


1.60 ± 0. 


,04 


1.10 


'b 


- Slogh 


< - 


-20 








4 


< 


z 


< 


0.6 


1.56 ± 0.16 


X 


10" 


-3 


12.50 ± 0.07 


0.37 ± 0. 


17 


10.88 ± 1.86 


13.51 ± 0.14 


1.18 ± 0.22 


1.68 ± 0. 


.04 


0.88 


' B 


- Slogh 


< - 


-20 


s 





4 


< 


" 


< 


0.6 


7.17 ± 0.78 


X 


10" 


-4 


12.77 ± 0.07 


0.30 ± 0. 


19 


4.97 ± 3.55 


13.88 ± 0.10 


1.18 ± 0.23 


1.83 ± 0. 


.06 


1.61 


'b 


- Slogh 


< - 


-19. 








.6 


< 


z 


< 


0.8 


3.88 ± 0.29 


X 


10" 


-3 


12.10 ± 0.07 


0.33 ± 0. 


19 


8.74 ± 5.33 


13.06 ± 0.12 


1.10 ± 0.15 


1.64 ± 0. 


,03 


0.58 


'b 


- Slogh 


< - 


-19. 


s 





6 


< 


z 


< 


0.8 


2.87 ± 0.22 


X 


10" 


-3 


12.23 ± 0.07 


0.37 ± 0. 


19 


12.47 ± 0.54 


12.96 ± 0.25 


0.82 ± 0.24 


1.69 ± 0. 


.04 


1.01 


' B 


- Slogh 


< - 


-20. 








6 


< 


z 


< 


0.8 


1.82 ± 0.15 


X 


10" 


-3 


12.38 ± 0.06 


0.26 ± 0. 


18 


5.79 ± 3.47 


13.40 ± 0.06 


1.12 ± 0.12 


1.77 ± 0. 


,04 


0.51 


' B 


- Slogh 


< - 


-20 


s 





6 


< 


z 


< 


0.8 


9.38 ± 0.85 


X 


10" 


-4 


12.64 ± 0.08 


0.34 ± 0. 


18 


3.96 ± 4.84 


13.69 ± 0.08 


1.16 ± 0.20 


1.91 ± 0. 


.05 


0.66 


' B 


- Slogh 


< - 


-21 








6 


< 


^ 


< 


0.8 


3.89 ± 0.37 


X 


10" 


-4 


12.93 ± 0.08 


0.28 ± 0. 


19 


10.01 ± 2.41 


14.07 ± 0.23 


1.33 ± 0.52 


2.15 ± 0. 


.08 


0.82 


' B 


- Slogh 


< - 


-19 


s 





8 


< 


z 


< 


1.0 


3.29 ± 0.24 


X 


10" 


-3 


12.09 ± 0.06 


0.24 ± 0. 


19 


11.21 ± 1.29 


13.13 ± 0.13 


1.22 ± 0.18 


1.80 ± 0. 


.03 


0.82 


' B 


- Slogh 


< - 


-20 








8 


< 


z 


< 


1.0 


2.26 ± 0.17 


X 


10" 


-3 


12.28 ± 0.07 


0.39 ± 0. 


17 


10.84 ± 1.90 


13.28 ± 0.12 


1.26 ± 0.18 


1.85 ± 0. 


.04 


1.09 


' B 


- Slogh 


< - 


-20 


s 





8 


< 


z 


< 


1.0 


1.22 ± 0.09 


X 


10" 


-3 


12.47 ± 0.07 


0.31 ± 0. 


18 


9.96 ± 1.49 


13.57 ± 0.07 


1.29 ± 0.17 


1.98 ± 0. 


.05 


1.12 


'b 


— Slogh 


< - 


-21 








8 


< 


^ 


< 


1.0 


5.34 ± 0.47 


X 


10" 


-4 


12.77 ± 0.07 


0.28 ± 0. 


19 


8.37 ± 3.29 


13.89 ± 0.12 


1.29 ± 0.32 


2.22 ± 0. 


.07 


1.28 


10 


-3-0fc3Mpc 


-3 







2 


< 


z 


< 


0.4 


1.00 ± 0.14 


X 


10" 


-3 


12.65 ± 0.09 


0.31 ± 0. 


19 


4.83 ± 3.88 


14.05 ± 0.25 


1.22 ± 0.56 


1.53 ± 0. 


.07 


1.31 


10 


-3-0fc3Mp^ 


-3 







4 


< 


z 


< 


0.6 


1.00 ± 0.12 


X 


10" 


-3 


12.64 ± 0.08 


0.31 ± 0. 


19 


12.36 ± 4.40 


13.67 ± 0.29 


1.00 ± 0.36 


1.75 ± 0. 


.06 


1.86 


10' 


-3-0).3mp^ 


-3 







6 


< 


z 


< 


0.8 


1.00 ± 0.09 


X 


10" 


-3 


12.59 ± 0.07 


0.25 ± 0. 


19 


10.41 ± 1.86 


13.68 ± 0.11 


1.10 ± 0.25 


1.90 ± 0. 


.05 


0.91 


10 


-SOh^Mpc 


-3 







8 


< 


^ 


< 


1.0 


1.00 ± 0.08 


X 


10" 


-3 


12.55 ± 0.07 


0.29 ± 0. 


19 


10.78 ± 2.25 


13.61 ± 0.10 


1.19 ± 0.27 


2.04 ± 0. 


.06 


1.17 


10' 


-^■'-•h^M-pc 


-3 







6 


< 


z 


< 


0.8 


3.16 ± 0.30 


X 


10" 


-4 


12.94 ± 0.07 


0.12 ± 0.16 


10.15 ± 0.80 


14.00 ± 0.15 


1.89 ± 0.35 


2.33 ± 0. 


.13 


0.87 


10' 


-3.5;,3[ 


vlpc 


-3 







8 


< 


z 


< 


1.0 


3.16 ± 0.28 


X 


10" 


-4 


12.93 ± 0.07 


0.27 ± 0.18 


9.06 ± 3.58 


14.20 ± 0.43 


1.16 ± 0.47 


2.38 ± 0.08 


0.63 



CO 



Table 5 
hod models of the luminosity and correlation functions of red galaxies with m^i^ and m[ as the only 



FREE PARAMETERS. 



sample Selection Photometric Space Density log Mj^^j^^ ctj^ j^ log Mq log m| a bg x /d.o.f- 

Criterion Redshift {h^ Mpc~^) 



'b 


- 51ogh < - 


-17.5 


0.2 


< 


z 


< 


0.4 


6.98 ± 0.92 


X 


10" 


-3 


11.94 ± 0. 


.05 


0.30 


11.94 ± 0.05 


12.73 ± 0.08 




.00 


1.44 ± 0. 


.04 


0.55 


'b 


- 51ogh < - 


-18.0 


0.2 


< 


z 


< 


0.4 


6.11 ± 0.80 


X 


10" 


-3 


11.98 ± 0. 


.05 


0.30 


11.98 ± 0.05 


12.81 ± 0.08 




.00 


1.43 ± 0. 


.04 


0.61 


'b 


- 51ogh < - 


-18.5 


0.2 


< 


z 


< 


0.4 


5.14 ± 0.69 


X 


10" 


-3 


12.03 ± 0.05 


0.30 


12.03 ± 0.05 


12.93 ± 0.09 




.00 


1.42 ± 0. 


.04 


0.66 


'b 


- 51ogh < - 


-19.0 


0.2 


< 


z 


< 


0.4 


3.96 ± 0.53 


X 


10" 


-3 


12.11 ± 0.05 


0.30 


12.11 ± 0.05 


13.10 ± 0.10 




.00 


1.40 ± 0. 


,04 


0.72 


'b 


- 51ogh < - 


-19.5 


0.2 


< 


z 


< 


0.4 


2.59 ± 0.36 


X 


10" 


-3 


12.27 ± 0.06 


0.30 


12.27 ± 0.06 


13.33 ± 0.11 




.00 


1.43 ± 0. 


,04 


0.75 


'b 


- Slogh < ' 


-20.0 


0.2 


< 


^ 


< 


0.4 


1.42 ± 0.20 


X 


10' 


-3 


12.50 ± 0. 


.06 


0.30 


12.50 ± 0.06 


13.72 ± 0.17 




.00 


1.48 ± 0. 


.05 


0.92 


'b 


- 51ogh < - 


-18.0 


0.4 


< 


^ 


< 


0.6 


5.75 ± 0.53 


X 


10"" 


-3 


12.00 ± 0. 


.04 


0.30 


12.00 ± 0.04 


12.74 ± 0.05 




.00 


1.57 ± 0. 


.02 


0.47 


'b 


- 51ogh < - 


-18.5 


0.4 


< 


z 


< 


0.6 


4.92 ± 0.46 


X 


10" 


-3 


12.05 ± 0. 


.03 


0.30 


12.05 ± 0.03 


12.82 ± 0.05 




.00 


1.58 ± 0. 


.02 


0.48 


'b 


- 51ogh < - 


-19.0 


0.4 


< 


z 


< 


0.6 


3.89 ± 0.37 


X 


10' 


-3 


12.14 ± 0. 


.04 


0.30 


12.14 ± 0.04 


12.93 ± 0.06 




.00 


1.60 ± 0. 


.02 


0.32 


'b 


- 51ogh < - 


-19.5 


0.4 


< 


z 


< 


0.6 


2.71 ± 0.27 


X 


10" 


-3 


12.26 ± 0. 


.04 


0.30 


12.26 ± 0.04 


13.14 ± 0.07 




.00 


1.61 ± 0. 


.03 


0.80 


'b 


- 51ogh < - 


-20.0 


0.4 


< 


z 


< 


0.6 


1.56 ± 0.16 


X 


10" 


-3 


12.47 ± 0. 


.04 


0.30 


12.47 ± 0.04 


13.41 ± 0.08 




.00 


1.69 ± 0. 


.03 


0.68 


'b 


- 51ogh < - 


-20.5 


0.4 


< 


^ 


< 


0.6 


7.17 ± 0.78 


X 


10" 


-4 


12.76 ± 0. 


.04 


0.30 


12.76 ± 0.04 


13.75 ± 0.09 




.00 


1.85 ± 0. 
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